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CHAPTER  1 

GENERAL  INFORMATION 


THE  FINITE  ELEMENT  CODE  ATILA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INST  I  TUT  SUPER  I  EUR  D* ELECTRON  I QUE  DU  NORD" 
(I  SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D' ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MAR I NE "  (GERDSM.  DCAN .  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE)  . 


Chapter  1.  General  Informat  Ion. 
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1.1  INTRODUCTION 


length 

bender 


The  finite  element  code  ATILA  has  been  specifically  developed 
to  aid  the  design  of  piezoelectric  devices,  mainly  for  sonar 
applications.  Thus,  It  Is  able  to  perform  the  modal  analysis  of 
both  ax Isymmetr ica I  and  fully  three-dimensional  piezoelectric 
transducers.  It  can  also  provide  their  harmonic  response  under 
radiating  conditions:  nearfleld  and  farfleld  pressure, 
transmitting  voltage  response,  directivity  pattern,  electrical 
impedance,  as  well  as  displacement  field,  nodal  plane  positions, 
stress  field  and  various  stress  criteria...  Its  accuracy  and  Its 
ability  to  describe  the  physical  behavior  of  various  transducers 
(Tonpllz  transducers,  double  headmass  symmetrical  length 
expanders,  free  flooded  rings,  f lextens Iona  I  transducers,  bender 
bars,  cylindrical  and  trilaminar  hydrophones...)  have  been  checked 
by  modelling  more  than  twenty  different  structures  and  comparing 
numerical  and  experimental  results.  These  tests  were  all 
successful,  provided  that  sufficient  ca^5f  was  previously  given  to 
the  physical  meanings  of  the  models.  Most\of  them,  as  well  as  the 
general  modelling  rules,  have  been  described  In  reports,  papers 
and  communications,  which  are  partial  I A  listed  In  the  next 
section.  \ 

ATILA  has  been  designed  for  the  "  Groups  d’ Etudes  et  de 
Recherche  en  Detection  Sous-Marlne,  G.E.R.D.S.M. ■  from  the 
"Direction  des  Constructions  et  Armes  Navales,  D.C.A.N.",  in 
TOULON  (FRANCE),  by  the  Acoustics  Laboratory  of  the  "Instltut 
Super ieur  d'Electronlque  du  Nord,  I.S.E.N.",  In  LILLE  (FRANCE).  It 
has  first  been  a  research  tool,  for  specific  optimisation 
purposes,  but,  due  to  the  request  for  Its  use  by  a  large  number  of 
people.  It  has  then  been  made  more  user  friendly  and  documented, 
while  a  specific  effort  has  been  carried  out  to  Improve  Its 
versatility  as  well  as  to  provide  more  efficient  graphic  display 
abilities.  At  the  present  time,  ATILA  has  been  run  on  various  IBM, 
VAX  and  UNIVAC  mainframes,  as  well  as  a  FPS  array  processor  and 
VAX  and  MASSCOMP  workstations.  Nevertheless,  ATILA  is  not  a  "black 
box"  and  users  are  assumed  to  have  a  firm  background  In  finite 
element  analysis. 

ATILA  has  been  developed  over  the  past  8  years  by  a  lot  of 
people,  but  mainly  by  K.  AN  I  FRAN  I,  R.  BOSSUT,  J.L.  CARTON,  J.C. 
DEBUS,  J.N.  DECARP  I GNY,  B.  DUBUS ,  R.  EDDE ,  B.  HAMONIC,  D.  MOREL, 
P.  TIERCE  at  I.S.E.N.,  and  B.  TOCQUET  and  D.  BOUCHER  at 
G.E.R.D.S.M..  Current  developments,  Including  new  algorithms,  new 
material  description,  material  loss  description  and  finite 
e iement/ l ntegra l  equations  coupling  are  proceeding  at  I.S.E.N., 
while  software  engineering  Is  now  provided  by  the  SINAPTEC 
company,  LILLE  (FRANCE). 


to  Improve  Its 
graphic  display 


^4^ 


Information  and  help  about  ATILA  can  be  obtained,  under 
specific  conditions,  from: 

.  J.N.  DECARP IGNY 
Acoustics  Laboratory 

Instltut  Superleur  d ' E lectron I que  du  Nord 
41  boulevard  Vauban 
59046  LILLE  Cedex .  FRANCE 

Tel:  (33) .20.30.62.20  Telex:  GEFIRN  120.369 

.  P.  TIERCE 
SINAPTEC  Company 
41  boulevard  Vauban 
59800  LILLE.  FRANCE 

Tel:  (33) .20.57.80.04  Telex:  GEFIRN  120.369 


Following  this  section  and  the  reference  list,  the  reader  will 
find  the  description  of  the  successive  stages  of  an  ATILA  Job  and 
their  connection  to  the  manual  organisation.  Then,  detailed 
Information  is  collected  in  five  chapters. 
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1.3  GENERAL  ORGANISATION  OP  AN  AT  I  LA  JOB 


Running  an  ATILA  Job  requires  different  well  defined  stages. 
This  section  describes  briefly  their  organisation  and  their 
connections  to  the  various  following  chapters  of  this  manual. 

1.3.1.  THE  MODEL  DEFINITION. 

Being  given  the  physical  problem,  the  user  has  to  define  the 
type  of  analysis  which  Is  required:  computation  of  the  static 
deformation  of  an  elastic  stucture  under  lumped  mechanical  load, 
modal  analysis  of  a  piezoelectric  stack  with  short  circuit 
electrical  condition,  computation  of  the  transmitting  voltage 
response  of  a  transducer...  This  type  of  analysis  has  to  be 
related  to  solving  algorithms  which  are  listed  and  described  in 
chapter  2.  New  solvers  are  currently  under  test  or  development  in 
l.S.E.N.  Acoustics  laboratory  and  detailed  Information  can  be 
provided  upon  request.  The  user  has  also  to  determine  the  type(s) 
of  elements  which  Is  (are)  required:  three-dimensional  elements, 
ax  I symmetr l ca I ,  plane  stress  or  plane  strain  elements,  plate  or 
shell  elements  relying  on  specific  approximations...  The  ATILA 
element  library  is  described  In  chapter  4. 

At  this  stage,  a  first  estimation  of  the  computation  size  can 
generally  be  made,  using  a  gross  determination  of  the  node  and 
degree  of  freedom  numbers  and  previous  CPU  time  measurements  or 
classical  numerical  rules.  This  estimation  Is  important  since  It 
can  avoid  the  generation  of  unrealistic  meshes.  At  this  stage 
also,  the  user  has  to  gather  the  requested  physical  parameters,  as 
described  In  chapter  4,  and  to  ensure  their  accuracy,  which  very 
definitely  determines  the  accuracy  of  the  result.  Though  not 
directly  Involved  In  the  data  file  generation,  this  stage  is 
fundamenta I . 

1.3.2.  THE  MESH  GENERATION. 

Following  the  previous  choice  of  element  type(s),  as  defined 
during  the  first  stage,  the  mesh  generation  implies  the  splitting 
of  the  whole  domain  Into  a  given  number  of  elements.  This  Is 
realised  by  defining  nodes  and  assigning  them  to  corresponding 
elements.  Node  definition  Is  made  by  listing  the  node  coordinates 
In  a  given  order,  called  the  node  numbering  order.  Then, 
assignment  of  the  nodes  to  the  corresponding  elements  is  obtained 
by  listing  for  each  element,  In  a  given  order  named  Its  topology, 
the  nodes  which  belong  to  It.  The  node  coordinate  description  Is 
detailed  In  chapter  3,  together  with  all  other  entries.  The 
topology  description  Is  presented  also  in  chapter  3,  after  the 
node  coordinate  description,  but  more  details  are  available  on  an 
element  by  element  basis  In  chapter  4. 
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During  the  mesh  generation,  the  user  can  use  the  Isoparametric 
character  of  all  the  ATILA  elements,  which  allows  easy  curved  line 
modelling  and  lower  density  meshes.  Nevertheless,  a  great  deal  of 
attention  has  to  be  paid  to  the  element  geometrical  aspect,  as 
well  as  to  the  mesh  step  which  has  to  be  related,  for  example,  to 
the  minimal  wavelength  value  of  Interest. 

Mesh  generation  can  be  fully  carried  out  by  the  user.  However, 
In  many  cases,  partial  regularity  of  the  mesh  allows  the  use  of  an 
automatic  generation  which  provides  node  coordinates  and  element 
topologies.  Then,  the  user  has  only  to  define  a  gross  splitting, 
built  with  super  elements,  and  to  select  their  automatic  splitting 
Into  finite  elements.  This  easier  procedure,  described  In  chapter 
6,  relies  upon  the  use  of  the  MOSAIQUE  code. 

1.3.3.  THE  DATA  FILE  PREPARATION. 

An  ATILA  data  file  Includes  entries,  specific  data  and 
comments.  It  Is  devoted  to  the  description  of  the  type  of 
analysis,  the  node  coordinates,  the  element  topology,  the  material 
properties,  the  element  geometrical  properties,  which  are  all 
entries,  and  the  loading  and  boundary  conditions  which  are 
formatted  data.  In  some  cases,  It  also  Includes  specific  mixed 
finite  element/plane  wave  model  data  and  graphic  display  data.  The 
whole  data  file  preparation  is  described  In  chapter  3  which  has  to 
be  carefully  read.  Specific  element  Information  which  are  required 
during  this  preparation  can  also  be  found  in  chapter  4.  Finally, 
when  the  automatic  mesh  generation  code  MOSAIQUE  Is  used,  data 
file  preparation  can  be,  In  a  large  part,  automatically  obtained, 
as  explained  In  chapter  6. 

A  detailed  check  of  the  data  file  can  be  processed  by  running  a 
specific  progam,  named  PGRAPH,  which  produces  a  complete  listing 
of  the  various  data.  Including  the  generated  degree  of  freedom 
numbering  which  reflects  exactly  the  prescribed  boundary 
conditions,  as  well  as  different  kinds  of  graphic  displays  (total 
or  partial  mesh,  with  or  without  element,  node  or  degree  of 
freedom  numberings).  The  use  of  PGRAPH  Is  described  in  chapter  5. 

1.3.4.  THE  RUNNING  OF  THE  CODE. 

Running  of  the  code  requires  the  edition  of  a  main  program 
which  has  to  set  the  array  sizes,  call  successive  subroutines  to 
read  the  data  file,  compute  the  elementary  stiffness  and  mass 
matrices,  assemble  these  matrices  into  the  global  stiffness  and 
mass  matrices,  solve  the  equations  and  display  the  results.  This 
progam  generation  Is  automatic,  given  the  corresponding  data  file, 
and  Is  provided  by  running  the  program  PGEN.  PGEN  will  Issue 
successive  questions,  on  the  terminal  screen,  to  specify  various 
user's  requests,  such  as  array  storages  for  a  delayed  solving, 
array  storages  for  a  postprocessing,  graphic  displays...  It  will 
provide  various  files,  Including  the  main  program  FORTRAN  source 
and,  In  the  case  of  a  VAX  computer  under  VMS  operating  system,  a 
COMMAND  file  to  compile  this  main  program  and  edit  the  links. 
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Running  of  the  code  will  provide  a  detailed  result  listing  and, 
according  to  the  specific  requests,  graphic  displays  or  stored 
arrays  for  restarting  another  solving  or  for  postprocessing.  It 
has  to  be  noted  that  PGEN  Is  also  able  to  generate  specific  main 
programs  to  restart  given  runs.  Every  Information  which  concerns 
the  program  running  as  well  as  PGEN  Is  contained  In  chapter  5. 

1.3.5.  RESULT  USING  AND  POSTPROCESSING. 

As  soon  as  the  run  Is  completed,  a  LISTING  file  Is  available 
which  contains  all  the  results:  elgenf requencles  and  eigenvectors, 
displacement  and  pressure  fields,  electrical  Impedance...  Tf 
requested,  specific  arrays  are  also  stored  for  postprocessing,  and 
mainly  the  array  PLO  which  contains  displacement  and  pressure 
fields  for  every  loading  case,  mode  number  or  excitation 
frequency.  Graphic  displays  of  the  distorted  structure  can  be 
easily  obtained  on  various  workstations  using  PGRAPH.  PGRAPH  will 
issue  FORTRAN  and  COMMAND  files  which  rely  upon  the  GKS  package 
and  have  to  be  executed  to  provide  the  display.  As  mentioned 
previously,  the  use  of  PGRAPH  Is  described  In  chapter  5. 

1 .3. 6. SUMMARY. 

As  a  summary,  the  next  page  displays  a  simple  flow  chart  of  an 
ATILA  finite  element  modelling.  A  large  part  of  this  flow  chart 
will  be  reproduced  In  a  more  detailed  and  accurate  fashion  In 
chapter  5,  with  reference  to  the  various  files  and  programs. 
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CHAPTER  2 
TYPES  OF  ANALYSIS 


THE  FINITE  ELEMENT  CODE  ATILA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  " INSTITUT  SUPERIEUR  D' ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D' ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MAR I NE "  (GERDSM.  DCAN .  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE)  . 


2.1  GENERAL  FORMULATION 


This  section  describes  the  general  formulation  which  has  been 
retained  to  model  radiating  piezoelectric  transducers  using  the 
ATILAoode  and  provides  a  complete  list  of  the  available  types  of 
analysis.  Then,  for  each  type  of  analysis,  the  following  sections 
give  more  Information  and  specific  details. 


Boundary 
surface  S. 


To  model  a  radiating  piezoelectric  transducer  using  the  ATILA 
code,  the  finite  element  mesh  must  Include  the  structure  as  well 
as  a  part  of  the  fluid  domain,  as  described  on  the  above  figure. 
Thus,  the  unknown  fields  are  the  displacement  field  u  In  the  whole 
structure,  the  electrical  potential  v  In  the  piezoelectric 
material  and  the  pressure  p  In  the  fluid.  Then,  the  equations 
which  are  solved  are,  first,  the  equation  of  motion  in  the 
structure,  second,  the  Poisson  equation  in  the  piezoelectric 
material  and,  third,  the  Helmholtz  equation  in  the  fluid. 
Kinematic  and  dynamic  continuity  equations  between  displacement 
and  pressure  fields  are  prescribed  at  the  Interface,  due  to  the 
variational  formulation,  and  an  appropriate  damping  condition  is 
applied  to  the  spherical  external  surface  S  that  surrounds  the 
fluid  doma I n . 


rvrrcwra _ 
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In  matrix  form,  the  complete  set  of  equations  associated  with 
this  problem  Is: 


[Kuu]  “  w2[M] 

EKUV] 

-  [L]  ;  !U!  IF 

III  1 

[Kuv]T 

m 

> 
> 
* 
i _ t 

co]T  i  iyi  -  !  -q 

III  1 

-a2c2w2[L]T 

[0] 

[H]  -  w2[Mi ] !  IP!  ! ac2F* 

where t 

-  U,  F,  V,  Q  and  P  are  vectors  that  contain  the  nodal  values  of 
the  displacement  field,  the  lumped  applied  forces,  the 
electrical  potential  and  charges,  the  pressure  field, 

-  F*  is  a  vector  that  contains  the  nodal  values  of  the  outgoing 
flux  through  S, 

-  [Kuu]  and  [H]  are  the  solid  and  fluid  stiffness  matrices, 

-  [M]  and  [Mi ]  are  the  solid  and  fluid  mass  matrices, 

-  [Kuv]  and  [Kvv]  result  from  the  assembling  of  the  piezoelectric 
and  dielectric  matrices, 

-  [L]  is  the  fluid-structure  interface  connectivity  matrix, 

-  [0]  Is  a  nul l  matr I x , 

-  a  and  c  are  the  fluid  density  and  sound  speed, 

-  w  is  the  circular  frequency, 

-  T  means  transposed. 

In  fact,  internal  electrical  potential  degrees  of  freedom  are 
associated  with  null  nodal  charges  and  can  be  easily  condensed. 
Moreover,  nodes  which  belong  to  hot  electrodes  have  the  same 
electrical  potential  degree  of  freedom  which  can  be  factored, 
following  simple  algebraic  operations.  Finally,  the  reference 
potential  can  be  arbitrarily  put  equal  to  zero.  Thus,  the 
preceedlng  matrix  equation  is  automatically  transformed  by  the 
ATILA  code  Into  the  following  form: 
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where  elastic,  piezoelectric  and  dielectric  terms  are  modified  bv 
the  Internal  electric  potential  static  condensation,  K  uv  and  0 
are  now  vectors,  Kvv  Is  now  a  scalar,  V  is  the  applied  electrical 
potential  and  I  is  the  electrical  current  entering  into  the 
structure.  Then,  If  a  dipolar  damping  condition  is  applied  on  the 
external  boundary  surface  S,  F*  Is  given  by: 


F*  -  -( 1/ac) . ( 1/R  +  Jw/c) .  [D]  P 

+  ( i/ac) . {( l/R  -  Jw/c)/(i  ♦  w2R2/c2)} . [O'  ]  P 
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Chapter  2.  Types  of  analysis 


where  R  Is  the  radius  of  S,  CD]  and  [D‘]  are  obtained  by 
assembling  radiating  surface  finite  elements  from  S,  the  first 
term  Is  the  monopolar  contribution  associated  with  a  spherical 
wave  Impedance  and  the  second  term  Is  the  dipolar  contribution. 

Starting  with  this  general  problem,  a  lot  of  different  specific 
analysis  can  be  processed.  Those  for  which  ATILA  provides  a  solver 
are  the  fol lowing: 

-  computation  of  the  displacement  field  of  an  elastic  structure 
under  a  static  lumped  or  distributed  load  (STA  1), 

-  computation  of  the  displacement  field  and  the  electrical 

potential  for  a  piezoelectric  structure  under  a  static  lumped  or 
distributed  load  (STA  2), 

-modal  analysis  of  an  elastic  structure,  providing  the 

elgenfrequencles  and  eigenmodes  (MOD  1), 

-modal  analysis  of  a  piezoelectric  structure,  providing  the 
elgenfrequencles  and  eigenmodes  under  resonance  or  ant l resonance 
electrical  conditions  (MOD  2), 

-  modal  analysis  of  a  closed  fluid  domain,  with  zero  pressure 

and/or  zero  flux  boundary  conditions,  providing  the 

elgenfrequencles  and  eigenmodes  (  MOD  3), 

-modal  analysis  of  a  closed  hydroelastic  system,  built  with  two 
coupled  elastic  and  closed  fluid  domains,  providing  the 
elgenfrequencles  and  eigenmodes  (MOD  4), 

-  modal  analysis  of  an  elastic  structure  flooded  in  an 

Incompressible  fluid  (added  mass  effect  only),  providing  the 
elgenfrequencles  and  eigenmodes  (MOD  5), 

-  harmonic  analysis  of  an  In-alr  piezoelectric  structure, 
providing  Its  electrical  Impedance  and  displacement  field  for  a 
given  excitation  frequency  (HAR  1), 

-  harmonic  analysis  of  an  In-alr  piezoelectric  structure,  using 
the  mixed  plane  wave/flnlte  element  method  and  providing  the 
displacement  field  for  a  given  excitation  frequency  (HAR  2), 

-  radiated  pressure  computation  for  a  given  displacement  field  of 
a  vibrating  structure  at  a  given  frequency  (HAR  3), 

-  harmonic  analysis  of  a  radiating  piezoelectric  transducer, 
providing  the  electrical  Impedance  and  the  displacement  and 
pressure  fields  In  the  whole  mesh  for  each  given  excitation 
frequency  (HAR  4), 

-  harmonic  analysis  of  a  radiating  piezoelectric  transducer,  using 
the  mixed  plane  wave/flnlte  element  method  and  providing  the 
displacement  and  pressure  fields  In  the  whole  mesh  for  each 
given  excitation  frequency  (HAR  5)- 


2.2  STATIC  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  (STA  1) 


In  this  case,  only  the  displacement  field  Is  concerned,  while  w 
has  to  be  put  equal  to  zero,  and  thus  the  whole  matrix  equation  Is 
reduced  to: 


ow  y  “  E 


The  user  has  to  provide  the  loading  F  and  the  solver  computes  U.  F 
has  to  be  defined  strictly  for  the”  finite  element  model,  which 
means  that,  for  example.  It  has  to  be  divided  by  2*pl  if  the  model 
is  ax  isymmetr  leal  ,  by  2  If  it  Is  applied  In  a  symmetry  plane 
limiting  the  mesh...  A  frontal  assembling  Is  performed,  coupled 
with  a  static  Gauss  condensation,  and  then  a  classical 
backsubst I  tut  Ion  algorithm  restores  the  whole  response.  Several 
different  loading  cases  can  be  solved  simultaneously.  The  user  has 
to  take  care  of  the  presence  of  rigid  body  modes,  which  are 
associated  with  singularities  of  the  system  of  equations.  Finally, 
If  dlstlbuted  loads  are  applied,  the  user  has  to  compute  the 


2.3  STATIC  ANALYSIS  OF  A  PIEZOELECTRIC  STRUCTURE  (STA  2) 


In  this  case,  the  matrix  equation  is  reduced  to: 
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If  the  external  electrical  connection  Is  a  short  circuit,  v  is 
equal  to  zero.  The  proceeding  equation  Is  reduced  again  to: 


ckuu]  y  -  f 


The  user  has  to  provide  the  loading  F  and  the  solver  computes  U, 
following  the  same  procedure  as  In  the  proceeding  section.  If  the 
external  electrical  connection  Is  open,  Q  Is  equal  to  zero  and: 
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The  user  has  to  provide  F  and  the  solver  computes  U  and  the 
Induced  electrical  potential  V.  Here  again,  attention  has  to  be 
paid  to  rigid  body  modes  as  well  as.  In  the  second  case,  to  the 
uniqueness  of  the  elect  I  cal  potential.  As  for  STA  1,  several 
loading  cases  can  be  processed  at  the  same  time. 
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2.4  MODAL  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  (MOD  1) 


In  this  case,  only  the  displacement  field  Is  of  concern  and  the 
loading  force  Is  put  equal  to  zero.  The  matrix  equation  Is  reduced 
to : 


aw  -  y  -  0 

Eigenvalues  and  eigenvectors  from  this  linear  system  of  equations 
are  the  elgenf requencles  and  elgenmodes  of  the  structure. 

The  ATILA  code  uses  frontal  assembling  coupled  to  a  Guyan 
condensation.  Thus,  the  user  has  to  provide  the  number  of  master 
degrees  of  freedom  which  have  to  be  retained  In  the  final 
d I agona I  I  sat  Ion  phase.  To  do  this,  he  has  to  keep  in  mind  that  the 
CPU  time  which  is  reqested  for  the  d I agona I  I  sat  ion  phase  Is 
roughly  proportional  to  the  cube  of  the  master  degree  number, 
absolute  value  of  this  time  depending  upon  the  algorithm,  the 
required  accuracy,  the  computer  and  the  compiler  option.  He  has 
also  to  keep  In  mind  that  accurate  e  I  genf  requenc  I  es  are  obtained 
only  for: 

w2  «  M I n(K | | /M | | ) 


where  K| (  and  M | (  are  the  diagonal  terms  of  the  stiffness  and  mass 
matrices  associated  with  the  condensed  slave  degrees  of  freedom. 
Thus,  he  has  to  perform  a  compromise  between  these  opposing 
requirements.  For  a  given  master  degree  number,  the  choice  of  the 
specific  degrees  which  have  to  be  retained  can  be  provided  by  the 
user  or  automatically  selected  by  the  code,  using  the  Henschell’s 
cr l ter i a . 

The  diagonal  I  sat  ion  algorithm  which  Is  used  by  ATILA  relies 
upon  the  classical  Householder  tr Idiagona I Isat Ion  method,  coupled 
to  a  Sturm  sequence.  It  can  be  processed  In  simple  or  double 
precision.  CPU  times,  as  obtained  on  an  IBM  4341  Mod.  2  (4  MO), 
are  reported  on  the  graphic  displays  of  pages  8  and  9.  Recently,  a 
new  d I agona I i sat  Ion  algorithm  has  been  Implemented.  It  provides 
the  e i genf requenc I es  in  double  precision  and  the  eigenvectors  in 
single  precision.  Though  much  more  efficient,  It  Is  still  under 
test  and  can  be  only  used  by  advanced  users,  following  a  short 
modification  of  the  ATILA  main  program*.  In  all  the  cases,  ATILA 
provides  all  the  el genf requencles,  the  number  of  which  is  equal  to 
the  number  of  master  degrees,  and  only  the  first  twenty 
elgenmodes,  if  available.  It  provides  also  the  rigid  body  modes, 
which  belong  to  the  previously  mentioned  first  twenty  modes.  Upon 
the  user's  request,  a  backsubst I  tut  Ion  routine  can  provide  the 
slave  degrees  of  freedom  for  any  of  the  obtained  elgenmodes. 
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In  some  cases,  care  has  to  be  taken  because  of  mixing  between 
modes  that  have  frequencies  which  are  too  close,  as  well  as  of 
matrix  I  I l-cond 1 1 Ion  I ng .  Successive  dlagona I Isat Ion  with  different 
precision  or  different  master  degree  numbers  can  Illuminate  these 
two  points. 


*•  To  use  the  new  d  I  agona  1  1  sat  Ion  algorithm,  the  user  has  to 
edit  the  main  program  FORTRAN  source,  as  obtained  for  a  double 
precision  eigenvalue  computation,  replace  CALL  D52200  (DLIST.S, 
S.CPDDC.PLO, IP.DM.Z)  by  CALL  S52200  (DL IST,S,CPDDC,PLO,  IP  ,Z)  and 
delete  the  DM  array  from  the  DIMENSION  set  of  instructions. 


on 


CPU  TIME  FOR  THE  ASSEMBLING  PROCESS  WITH  RESPECT  TO  THE  NUMBER  OF 
MASTER  DEGREES  OF  FREEDOM.  The  parameter  is  the  total  number  of 

degrees  of  freedom. 
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2.5  MODAL  ANALYSIS  OF  A  PIEZOELECTRIC  STRUCTURE  (MOD  2) 


In  this  case,  the  matrix  equation  Is  reduced  to: 


!  [Kuu]  "  *2[M] 
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If  a  resonance  condition  occurs,  V  Is  equal  to  zero  and: 


{[Kuu]  -  w2[M]}  U  -  0 

If  an  ant  I  resonance  condition  occurs,  I 
{[Kuu]  -  w2CM]  -  K  UV.K  uvt/Kvv}  U  - 


Is  equal  to  zero  and: 

0 


In  the  two  cases,  all  the  Information  given  In  the  proceed ing 
section  with  respect  to  the  modal  analysis  of  an  elastic  structure 
is  applicable.  The  only  point  which  has  to  be  added  is  that,  due 
to  the  automatic  static  condensation  of  the  electrical  potential 
degrees  of  freedom,  all  these  degrees  of  freedom  are  slaves. 


2.6  MODAL  ANALYSIS  OF  A  CLOSED  FLUID  DOMAIN  (MOD  3) 


In  this  case,  only  the  pressure  field  Is  concerned.  Moreover, 
as  zero  pressure  and/or  zero  flux  boundary  conditions  have  to  be 
applied,  the  F'  components  which  are  able  to  appear  In  the 
equation  are  all  zero.  Thus: 


{[H]  -  w2^]}  P  -  0 

All  the  Information  provided  In  section  2.5  with  respect  to  the 
modal  analysis  of  an  elastic  structure  Is  applicable.  Moreover, 
the  user  has  to  remember  that  the  zero  flux  condition  Is  the 
natural  boundary  condition  of  this  finite  element  formulation  and 
automatically  applies  to  any  boundary,  If  not  superseded  by 
another  one.  This  modal  analysis  can  be  used  to  determine  the 
Irregular  frequencies  associated  with  an  exterior  Helmholtz 
Integral  equation. 


2.7  MODAL  ANALYSIS  OF  AN  HYDROELASTIC  SYSTEM  (MOD  4) 


In  this  case,  structure  and  fluid  domains  exist  together  and 
Interact  along  their  Interface.  However,  the  fluid  domain  has  to 
be  closed,  no  radiation  damping  being  Included  In  the  model.  Thus, 
fluid  surfaces  which  are  not  In  contact  with  the  solid  structure 
have  to  be  zero  pressure  or  zero  flux  surfaces,  as  Is  the  case, 
for  example,  for  fluid  In  an  open  tank.  Restricting  the  general 
matrix  equation  of  section  2.1  to  the  displacement  and  pressure 
field,  the  equation  which  Is  obtained  Is: 


[kuu]  -  w2tMl  -  ^ 


U 


0 


-a2c2w2[l]T 


[H]  -  w2[M1] 


P 
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Apart  for  a  slight  difference  regarding  the  algorithm,  all  the 
Information  given  In  section  2.4  with  respect  to  the  modal 
analysis  of  an  elastic  structure  Is  applicable.  With  respect  to 
the  algorithm,  the  above  matrix  has  to  be  made  symmetrical  before 
the  ATILA  classical  d I agona I  I  sat  Ion  is  performed.  This 
symmetr Isat Ion  is  processed  using  a  standard  method  due  to  Irons. 
However,  It  has  to  be  pointed  out  that  the  symmetr  i  sat  Ion  method 
requires  the  Inversion  of  [H]  and  then  Implies  that  no  fluid 
Isopressure  mode,  equivalent  to  structure  rigid  body  modes, 
exists.  Thus,  at  least  a  small  part  of  the  fluid  surface  has  to  be 
a  zero  pressure  surface.  Recent  developments  have  led  to  the 
design  of  a  new  solver  which  avoids  this  difficulty^  but  it  Is 
still  under  test  and  can  only  be  run  by  advanced  users  .  Moreover, 
during  the  assembling  phase,  as  well  as  for  an  elastic  structure, 
Guyan  condensation  Is  performed  for  both  solid  and  fluid  degrees 
of  freedom.  Care  has  to  be  given  to  the  Interface  degrees  of 
freedom  which  have  to  be  kept  in  sufficient  number.  Finally, 
piezoelectric  parts  can  be  Included  In  the  structure.  Then,  the 
user  has  also  to  provide  the  corresponding  electrical  boundary 
cond l t Ions . 


*.  To  use  this  new  algorithm,  the  user  has  to  edit  the  main 
program  FORTRAN  source,  as  obtained  for  this  analysis,  delete  CALL 
R52300  and  replace  CALL  R52200  by  CALL  S52300  with  the  same 
parameters.  The  new  algorithm  detects  automatically  the 
Isopressure  or  rigid  body  modes  and  deals  with  them  using  a 
deflation  method.  In  opposition  to  the  standard  algorithm,  It  can 
also  compute  the  elgenmodes  for  mere  fluid  or  solid  domains, 
without  Interaction. 
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2.8  MODAL  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  FLOODED  IN  AN 
INCOMPRESSIBLE  FLUID  (MOD  5) 


Incompressible  fluid  condition  is  an  approximation  that  Is 
useful  to  compute  the  fluid  mass  loading  effect  on  the  low 
frequency  modes  of  an  elastic  structure.  Obviously,  It  discards 
any  radiation  effect.  If  It  Is  applied  to  the  matrix  equation  of 
the  previous  section,  this  equation  becomes: 


C(KUU]  -  w2([M]  +  a2c2[L][H]-1[L]T)}  U  -  0 


where  a2c2[L] [H]-1 [L]T  is  the  added  mass.  Then,  the  classical 
ATILA  dlagona I isat Ion  provides  the  structure  elgenfrequencies  and 
eigenmodes.  To  obtain  this  particular  condition  using  ATILA 
requires  that  all  the  fluid  degrees  of  freedom  are  condensed,  this 
approximation  producing  exactly  the  same  result  as  the 
incompressibility  hypothesis*.  If  the  fluid  domain  extension  Is 
Infinite,  It  has  to  be  simply  truncated  before  the  mesh 
generation.  Successive  truncations  at  various  distances  from  the 
wet  structure  are  the  best  convergence  test. 


*.  After  the  Guyan  condensation  of  all  the  fluid  degrees  of 
freedom,  the  matrix  which  has  to  be  dlagonailsed  Is  already 
symmetrical  but  the  automatic  program  generator  will  provide  the 
call  for  the  symmetr I  sat  Ion  subroutine  because  the  data  file 
corresponds  to  a  fluid-structure  modal  analysis.  Thus,  before 
running  ATILA  the  user  must  edit  the  FORTRAN  file  and  delete  the 
I  ine  CALL  R52300 . 
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2.9  HARMONIC  ANALYSIS  OF  AN  IN-AIR  PIEZOELECTRIC  STRUCTURE  (HAR  1) 


In  this  case,  an  electrical  potential  of  given  frequency  Is 
applied  to  an  In-alr  free  piezoelectric  structure.  Then,  the 
general  matrix  equation  of  section  2.1  can  be  reduced  to: 


ICKuul  -  w2£m3 
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To  obtain  this  whole  equation  using  the  ATILA  code  Implies  that 
the  applied  electrical  potential  degree  of  freedom  Is  not 
automatically  condensed,  which  requires  a  specific  Identification 
using  the  boundary  condition  data.  Then,  the  solver  computes  the 
displacement  field  U,  using  the  first  line  of  the  above  matrix 
equat ion: 


{[Kuu]  ~  w2[M]}  U  -  -K  uv.V 


as  well  as  the  electrical  impedance  Z,  using  the  second  line: 


1/Z  -  -Jw  (Kvv  +  K  uvT.U/V) 


The  solving  algorithm  Is  a  simple  Gauss  elimination.  In  the  same 
run,  solving  can  be  performed  for  20  different  excitation 
frequencies.  Moreover,  Job  restart  Is  possible  for  new  frequency 
values  without  performing  again  the  whole  assembling.  The 
Impedance  value  which  Is  provided  corresponds  strictly  to  the 
finite  element  model.  Thus,  the  user  has  to  divided  the  result  by 
2*pl  If  the  model  is  ax  I  symmetr  I  ca  I  ,  2  If  the  mesh  Is  restricted 
to  half  of  the  structure  by  a  symmetry  plane,  etc. 
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2.10  HARMONIC  ANALYSIS  OF  AN  IN-AIR  PIEZOELECTRIC  STRUCTURE  USING 

THE  MIXED  METHOD  (HAR  2) 


In  this  case,  the  driving  part  of  the  structure  Is  modelled 
using  the  plane  wave  approximation  and  a  transfer  matrix  method, 
while  the  remaining  Inactive  part  of  this  structure  Is  modelled 
using  finite  element.  These  two  parts  are  In  contact  along  a  plane 
Interface,  which  exhibits  also  the  plane  wave  type  displacement 
field.  Thus,  all  the  degrees  of  freedom  which  belong  to  the  finite 
element  side  of  the  interface  and  correspond  to  a  mechanical 
displacement  normal  to  it  have  to  be  Identical  and  are  denoted 
hereafter  as  UQ.  Then,  the  transfer  matrix  method  provides,  for 
each  given  frequency,  a  mechanical  force  F0  which  Is  generated  by 
the  piezoelectric  effect  and  Is  applied  to  the  finite  element  part 
at  the  Interface.  This  force  can  be  written: 


Fq  “  Zq  .  Uq  +  Fy 


where  ZQ  Is  a  mechanical  Impedance  and  Fv  contains  the  electrical 
excitation  effect.  Finally,  the  general  matrix  equation  becomes: 


,TKUU]  -  w2[M] 
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where  [Kuu]  and  [M]  correspond  to  all  the  degrees  of  freedom  apart 
from  UQ .  The  user  has  to  select  the  excitation  frequencies  and  to 
specify,  using  the  correct  boundary  condition  data,  the  degree  of 
freedom  UQ.  The  solving  algorithm  Is  a  simple  Gauss  elimination, 
and  the  code  provides  the  displacement  field  In  the  finite  element 
part  of  the  structure.  Twenty  different  excitation  frequencies  can 
be  provided  for  the  same  run.  Moreover,  as  In  section  2.9,  Job 
restart  Is  possible  with  new  frequency  values.  This  mixed  method 
Is  generally  well  suited  to  the  modelling  of  Tonpilz  type 
transducers  and  allows  large  savings  of  CPU  time  . 


2.11  RADIATED  PRESSURE  COMPUTATION  FOR  A  GIVEN  DISPLACEMENT  FIELD 

(HAR  3) 


In  this  case,  a  vibrating  structure  radiates  an  acoustic  wave 
in  an  Infinite  medium  and  the  vibration  displacement  Is  assumed  to 
be  known  for  each  given  frequency.  The  general  matrix  equation 
from  section  2.1  can  be  reduced  to: 

C [H]  -  W2[M1]}  P  -  ac2F '  +  a2c2w2[L]T.U 


Then,  If  monopolar  or  dipolar  damping  Is  used,  following  the  third 
equation  of  section  2.1,  the  damping  term  can  be  written  under  the 
form : 


ac2F'  -  [G]  P 


where  [G]  stands  for  a  linear  complex  frequency  dependent 
operator.  Thus: 

CCH]  -  w2^]  -  [G])  P  -  a2c2w2[L]T.u 


The  user  has  to  define  U  and  the  given  excitation  frequencies. 
To  define  U,  he  has  to  write  a  small  subroutine,  named  FABU* , 
which  contains  the  displacement  value  for  each  Interface  solid 
degree  of  freedom,  and  has  to  be  added  to  the  main  program.  Only 
one  displacement  field  can  be  defined  at  a  time  and  this  field  has 
to  be  real  (phase  reference).  Ten  successive  frequencies  can  be 
processed  In  the  same  run  and,  as  for  the  previous  section, 
restart  is  available.  The  solver  splits  the  previous  complex 
equation  Into  two  real  equations.  Then,  a  matrix  Inversion 
followed  by  a  Gauss  elimination  is  performed.  For  each  frequency, 
ATILA  provides  the  real  and  Imaginary  parts  of  the  pressure  field, 
printed  in  two  successive  tables. 

To  obtain  accurate  results,  the  user  has  to  give  a  great  deal 
of  attention  to  the  damping  element  choice.  If  monopolar  damping 
Is  used,  the  external  boundary  surface  has  to  be  centered  at  the 
acoustic  center  and  it  has  to  be  In  the  farfleld.  If  dipolar 
damping  Is  used,  the  external  boundary  surface  has  also  to  be 
centered  at  the  acoustic  center  but,  generally,  It  can  be  deeply 
In  the  nearfleld.  However,  In  this  second  case,  farfleld 
Information  Is  only  available  after  an  appropriate  postprocessing. 
Dipolar  elements  can  only  be  used  for  ax Isymmetr  lea  I  models. 
Three-dimensional  dipolar  element  design  Is  In  progress. 
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Degrees  of  freedom  which  belong  to  the  exterior  boundary 
surface  are  automatically  master  degrees.  Those  which  belong  to 
the  Interface  can  be  slave  degrees  but  the  user  has  to  keep  a 
sufficient  number  of  them  as  master  degrees.  Finally,  care  has  to 
be  taken  to  keep  the  maximal  distance  between  two  successive  fluid 
master  nodes  In  the  propagation  direction  less  than  one  quarter  of 
a  wavelength. 


* .  The  FORTRAN  subroutine  which  has  to  be  added  at  the  end  of 
the  main  program,  to  provide  the  given  displacement  field  U,  Is 
FABU  (CPDDC.U, IOLE, IP) .  It  Is  called  by  the  solving  subroutine 
R52400  during  the  run.  IOLE  and  IP  are  control  Integers  which  are 
defined  In  R52400,  as  well  as  the  CPDDC  array.  In  return,  FABU 
provides  all  the  components  of  U,  numbered  In  the  ascending  order. 
U  has  to  contain  IDT  components,  where  IDT  is  the  total  number  of 
degrees  of  freedom,  but,  In  fact,  only  the  mechanical  degrees  of 
freedom  are  concerned.  The  following  lines  are  provided  as  a 
standard  example: 

SUBROUTINE  FABU ( CPDDC ,U , I  OLE , I P ) 

DIMENSION  CPDDCdP,*),  U(») 

U(1)  -  ... 

U( 2)  -  ... 


RETURN 
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2.12  HARMONIC  ANALYSIS  OF  A  PIEZOELECTRIC  RADIATING  STRUCTURE 

(HAR  4) 


This  Is  the  general  case.  Using  the  operator  [G]  defined  In  the 
previous  section  2.11,  the  matrix  equation  Is: 
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and: 


1/Z  -  -Jw  (K  UVT.U  +  Kvv . V) 


The  user  provides  the  excitation  frequencies  and,  using 
specific  boundary  condition  data,  defines  the  degree  of  freedom 
that  Is  associated  with  the  applied  electrical  potential.  The 
solver  computes,  using  matrix  Inversion  and  Gauss  condensation, 
the  real  and  Imaginary  parts  of  the  displacement  and  pressure 
fields.  Then,  it  computes  the  electrical  impedance.  Real  and 
Imaginary  parts  are  provided  in  two  successive  tables.  Ten 
successive  excitation  frequencies  can  be  processed  In  the  same 
run,  and,  as  In  the  previous  section,  restart  Is  available. 

Information  given  In  the  previous  section  and  regarding  damping 
elements,  nearfleld  or  farfield  conditions,  boundary  surface  and 
interface  degrees  of  freedom  remain  available  for  this  section  and 
are  not  duplicated.  Information  given  In  section  2.9  regarding  the 
electrical  Impedance  Is  also  valid. 


2.13  HARMONIC  ANALYSIS  OF  A  PIEZOELECTRIC  RADIATING  STRUCTURE 
USING  THE  MIXED  METHOD  (HAR  5) 


This  case  Is  a  mixing  of  the  previous  cases  reported  In 
sections  2.10  and  2.12.  Using  the  corresponding  notations,  the 
matrix  equation  becomes: 
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The  user  provides  the  excitation  frequencies  and,  using 
specific  boundary  condition  data,  defines  the  degree  of  freedom 
that  is  associated  with  the  normal  displacement  to  the  Interface. 
The  solver  computes,  using  matrix  Inversion  and  Gauss 
condensation,  the  real  and  Imaginary  parts  of  the  displacement  and 
pressure  fields.  Real  and  Imaginary  parts  are  provided  in  two 
successive  tables.  Ten  successive  excitation  frequencies  can  be 
processed  In  the  same  run,  and,  as  In  the  previous  section, 
restart  Is  aval  I  able. 

Information  given  In  section  2.11  regarding  damping  elements, 
nearfleld  or  farfleld  conditions,  boundary  surface  and  Interface 
degrees  of  freedom  remain  available  for  this  section  and  are  not 
dupl icated. 
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THE  FINITE  ELEMENT  CODE  ATILA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INSTITUT  SUPERIEUR  D' ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D' ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN .  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE)  . 


3.1  INTRODUCTION 


This  section  Is  placed  first  In  the  data  preparation  part  of 
this  manual  for  pedagogical  purposes  only.  In  fact,  In  most  cases, 
an  automatic  data  generation  can  be  processed  using  the  MOSAIQUE 
code  and  a  more  simple  Initial  data  file.  This  automatic  data 
generation  Is  described  in  the  chapter  6.  When  used  first, 
MOSAIQUE  provides  data  files  which  are  Identical  to  the  files 
described  hereafter. 

in  the  general  case,  an  ATILA  data  file  Includes  five  different 
sect  ions : 

-  the  header,  which  contains  only  comments  and  describes  the  Job 
under  consideration, 

-  a  set  of  entries,  which  allows,  among  other  things,  to  describe 
the  mesh  and  the  properties  of  the  materials, 

-  a  set  of  data,  which  describes  the  loading  and  boundary 
cond I t Ions , 

-  a  set  of  data,  which  describes  the  plane  wave  part  of  the 
structure  In  the  case  of  a  mixed  model, 

-  a  set  of  data,  which  defines  the  graphic  displays. 

All  data,  except  for  the  angles,  must  be  given  using  MKS  units  or 
any  coherent  unit  system  which  Is  deduced  from  the  MKS  system.  The 
angle  values  are  always  given  In  degrees. 

In  this  chapter,  we  always  refer  to  a  simple  Illustrative 
example  and  give  the  corresponding  version  of  each  entry  or  data 
set.  This  example  is  the  modal  analysis  of  an  elementary 
ax i symmetr l ca l  piezoelectric  structure,  made  up  of  two  ceramic 
rings  and  two  thick  electrodes,  under  electrical  ant  I  resonance 
conditions.  This  structure  has  a  symmetry  plane  which  Is  normal  to 
the  axis  and  the  modes  which  have  to  be  computed  are  symmetrical 
with  respect  to  this  plane.  Five  full  eigenvectors  are  requested. 
The  structure,  its  material  properties  and  the  corresponding  mesh 
are  described  on  the  following  two  pages.  Other  examples  are  given 
at  the  end  of  the  manual. 


1 


% 
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4esh. 


.  1 


2  .  . 


SYMMETRY 

AXIS 


Due  to  simple  symmetry  arguments,  the  mesh  Is  limited  to  one  half 
of  the  structure.  Thus,  the  OY  axis  Is  the  Intersection  of  the 
symmetry  plane  with  the  plane  of  the  figure. 


Chapter  3.  Description  of  a  data  file. 


3.2  THE  HEADER 


The  header  includes  comments  which  describe  the  job  under 
consideration.  The  whole  header  Is  automatically  copied  on  the 
first  page  of  every  listing  which  is  related  to  the  job.  The  first 
line  of  the  header,  generally  a  title,  Is  automatically  repeated 
ahead  of  each  listing  page.  To  define  the  header,  each  line  begins 
with  the  symbol  *.  The  number  of  lines  In  a  header  Is  unlimited. 


EXAMPLE. 


+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


* 

* 

* 


.  MODAL  ANALYSIS  .  PIEZOELECTRIC  STRUCTURE  . 


.  ANT  I  RESONANCE  CASE  . 

.  AX  I SYMMETR I  CAL  MODEL  . 

.  1/2  OF  THE  STRUCTURE  (SYMMETRY  PLANE)  _ 


.  NUMBER  OF  ELEMENTS:  2  .. 

.  NUMBER  OF  NODES:  13  _ 

.  NUMBER  OF  DOF:  26  . 

.  NUMBER  OF  MASTER  DOF:  23 
.  DOUBLE  PRECISION  . 


.  DATA  FILE  WRITTEN  BY: 
.  DATE  . 


3.3  THE  ENTRIES 


The  various  entries  which  are  described  hereafter  are  listed  in 
a  given  order,  which  corresponds  to  a  classical  use  but  Is  not 
mandatory.  These  entries  have  to  be  written  using  a  simple  free 
format,  described  below,  which  allows  an  easy  reading  of  the  file. 
An  entry  can  be  associated  to  given  peculiar  words,  named 
attributes,  which  have  to  be  written  on  the  entry  line,  or  to 
given  numerical  values,  named  parameters,  which  have  to  be  written 
on  the  following  line.  Entries  being  read  In  free  format,  all  the 
attributes  and  parameters  have  to  be  provided,  even  If  they  are 
dummy  or  if  their  numerical  value  Is  zero. 


3.3.1  DESCRIPTION  OF  THE  FREE  FORMAT. 

Control  characters. 

.  [  ]  and  [,]  separate  two  words  on  the  same  line.  They  have  to  be 
used  to  separate  an  entry  from  Its  following  attribute,  to 
separate  two  following  attributes  or  to  separate  two  following 
parameters . 

.  [/]  and  [-]  separate  two  following  lines  (carriage  return).  As 
well  as  the  carriage  return,  they  have  to  be  used  to  separate  an 
entry  from'lts  following  parameter. 

.  [7]  deletes  the  preceedlng  characters  of  the  line, 

.  [»]  when  typed  at  the  beginning  of  a  new  line,  means  that  this 
line  Is  a  comment.  This  line  is  not  printed  on  a  listing.  Such  a 
comment  has  not  to  be  confused  with  the  header. 

.  [* . *]  allows  the  writing  of  a  comment  in  an  entry 

line.  This  comment  is  not  printed  on  a  listing. 

Real  number  syntax. 

The  following  syntaxes,  given  as  examples,  are  available: 

1.22E+02  .  -1.22E2 

-IE-3  .  1.22 

-1  .  -1 

Examp  I es . 

The  following  sequences  are  equivalent  to  define  the  given 
example:  one  alphanumeric  datum  (alpha)  on  the  first  line,  two 
Integers  (1,2)  on  the  second  line,  three  real  numbers  (1..2..3.) 
on  the  third  line: 

.  alpha  .  alpha-1,2  .  alpha/l  2-1..2.  3. 

1,2  1 .  2. ,3. 

1 . ,2. ,3. 


3.3.2  THE  ENTRY  DIRECTORY 


NAME 

COMP. 

ATTRIBUTES  OR  PARAMETERS 

P 

ANALYSIS 

YES 

STATIC,  MODAL  or  HARMONIC 

7 

CLASS 

AX  1 SYMMETR 1  CAL ,  PSTRESS  or 
PSTRAIN 

1 

NLOAD 

YES 

NLO 

7 

PRECISION 

NO 

SINGLE  or  DOUBLE 

8 

LCPDDC 

YES 

NLC 

8 

NODES 

YES 

node  coordinates 

(end  of  the  set:  dummy  line) 

9 

NEWAXES 

NO 

CARTESIAN,  CYLINDRICAL  or 
SPHERICAL 

XO  YO  ZO  OZZ  OYY  OXX 

10 

SCALE 

NO 

EX  EY  EZ 

11 

ELEMENTS 

■ 

Rxxxxx  MATERx  NGx 

Element  topology 

(end  of  the  set:  dummy  line) 

12 

MATERIALS 

YES 

MATERx 

Properties  of  the  material. 

(end  of  the  set:  dummy  line) 

14 

GEOMETRY 

NO 

NGx 

Geometrical  properties. 

(end  of  the  set:  dummy  line) 

17 

GEOMETRY 
POLAR IZAT. 

NO 

CARTESIAN,  CYLINDRICAL  or 
SPHERICAL. 

NGx 

Euler  angles,  coordinates. 

(end  of  the  set:  dummy  line) 

18 

REDUCTIONS 

YES 

1  N  MDF 

24 

FREQUENCIES 

NO 

Frequency  values. 

25 

RADIATION 

NO 

MONOPOLAR  or  DIPOLAR. 

25 

PRINTING 

NO 

NPR 

25 

END 

YES 

26 

COMP,  means  that  the  entry  Is  compulsory. 


fi 


3.3.3  DESCRIPTION  OF  THE  ENTRIES. 


In  this  section,  the  use  of  the  symbol  9  In  front  of  the  entry 
name  means  that  this  entry  Is  compulsory  and  has  to  be  used  In 
every  Job.  The  other  entries  have  to  be  used  only  In  specific  Jobs 
or  are  opt  Iona  I . 


«  ANALYSIS  STATIC  or  MODAL  or  HARMONIC 

This  entry  defines  the  kind  of  analysis  which  Is  requested  in  the 
job  (see  chapter  2,  theoretical  background)  and  Is  compulsory. 


Test  example.  + 

+  ANALYS I S  MODAL 
+ 


CLASS  AX  I  SYMMETRICAL  or  PLSTRESS  or  PLSTRAIN 

This  entry  specifies.  If  necessary,  that  the  analysis  corresponds 
to  an  ax isymmetr lea l  ,  a  plane  stress  or  a  plane  strain  model.  It 
is  only  needed  If  2D  elastic  elements  are  used,  but  then  Is 
compulsory. 


Test  example.  + 

+  CLASS  AX  I  SYMMETRICAL  . 


•  NLOAD 
NLO 

This  entry  Is  compulsory  and  provides  NLO  which  Is: 

.  the  number  of  loading  cases  for  a  static  analysis, 

.  the  number  of  full  eigenvectors  which  are  required  In  the  case 
of  a  modal  analysis, 

.  the  number  of  given  frequency  values  In  the  case  of  an  in-alr 
harmonic  analysis, 

.  twice  the  number  of  given  frequency  values  In  the  case  of  an 
harmonic  analysis  of  a  radiating  structure. 


Test  example. 


+ 

+  NLOAD  -  5 


PRECISION  SINGLE  or  DOUBLE 


This  entry  is  only  used  If  an  eigenvalue  computation  is  run  after 
the  assembling  and  specifies  the  precision  of  this  computation. 
The  default  answer  is  single. 


Test  example.  + 

+  PRECISION  DOUBLE 
+ 


a  LCPDDC 
NLC 


The  CPDDC  array  contains  the  node  coordinates  (X,  Y,  Z)  and  the 
numbers  of  the  node  degrees  of  freedom  (translations  UX,  UY,  UZ , 
rotations  NX,  NY,  NZ,  pressure  P  and  electrical  potential  V) . 
These  data  are  listed,  with  one  row  per  node.  In  the  following 
order : 

X.  Y,  Z,  UX  (or  P),  UY,  UZ,  NX  (or  V) .  NY,  NZ 

The  LCPDDC  entry  is  compulsory  and  provides  NLC,  which  is  the 
number  of  terms  needed  in  one  row.  Thus,  In  the  case  of  an 
ax Isymmetr lea  l  analysis  of  an  elastic  structure  NLC  equals  5, 
while  In  an  ax  I symmetr I ca I  analysis  of  a  piezoelectric  structure 
NLC  equals  7  and  In  an  analysis 'of  a  plate  NLC  equals  8.  If 
useless  Intermediate  columns  are  Introduced,  due  to  the  given 
value  of  NLC,  the  corresponding  degrees  of  freedom  will  have  to 
be  deleted  later.  This  case  happens,  as  an  example,  for  the 
degrees  of  freedom  Included  In  columns  3,  4  and  5  in  the  plate 
case. 


+ 

+  LCPDDC 
+  7  .  .  .  . 
+ 


Test  example. 
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XI  Y1  Z1 
X2  Y2  Z2 
X3  Y3  Z3 


XN  YN  ZN 
.DUMMY. . 


This  entry  provides  the  node  coordinates,  listed  with  one  row  per 
node.  In  the  node  number  ascending  order.  The  list  has  to  be  ended 
by  a  dummy  line  and  the  entry  Is  compulsory.  If,  for  any  reason, 
the  node  coordinate  list  Is  spilt  in  several  successive  sets,  each 
set  requires  a  new  entry  NODES  and  has  to  be  finished  by  a  new 
dummy  line. 


_  WARNING  _ 

When  the  model  Includes  both  structure 
and  fluid  domains,  nodes  have  to  be  set 
In  the  following  order: 

SOLID  -  FLUID  -  RADIATING. 

For  ax Isymmetr leal  models,  the  OX  axis 
Is  necessarily  the  symmetry  axis. 


Test  example. 


+  NODES 
+  0.00 
+  0.00 
+  0.00 
+  0.01 
+  0.01 
+  0.02 
+  0.02 
+  0.02 
+  0.04 
+  0.04 
+  0.06 
+  0.06 
+  0.06 


NEWAXES  CARTESIAN  or  CYLINDRICAL  or  SPHERICAL 
XO  YO  ZO  OZZ  OYY  OXX 

The  Initial  coordinate  system  In  which  the  node  coordinates  are 
given  (default  system)  is  cartesian  and  Its  origin  and  axes  are 
the  absolute  origin  and  axes.  The  entry  NEWAXES  allows  it  to  be 
mod  If  led. 

If  the  new  system  Is  cylindrical  (spherical),  the  node  coordinates 
X,  Y,  Z,  associated  to  the  following  NODES  entry(les)  become  the 
cylindrical  (spherical)  coordinates  R,  THETA,  Z  (R,  THETA,  PHI). 
Moreover,  XO,  YO,  ZO  are  the  coordinates  of  the  origin  of  the  new 
system,  while  OZZ,  OYY,  OXX  are  the  angle  values  of  the  Euler 
rotations  which  transform  the  Initial  axes  In  the  new  axes.  The 
first  rotation  OZZ  Is  around  the  OZ  axis  of  the  Initial  system. 
The  second  rotation  OYY  is  around  the  OY  axis  provided  by  the 
first  rotation.  The  third  rotation  OXX  Is  around  the  OX  axis 
provided  by  the  second  rotation.  Angle  values  have  to  be  given  in 
degrees . 

This  entry  can  be  used  several  times  during  the  description  of  the 
mesh,  Its  parameters  being  always  defined  with  respect  to  the 
absolute  system.  To  come  back  to  the  absolute  system,  It  has  to  be 
used  with  the  CARTESIAN  attribute,  the  parameters  being  equal  to 
zero. 


Test  example.  + 

+  NODES  .... 
+0.00  0.02 
+0.00  0.05 
+0.00  0.08 
+  0.01  0.02 
+0.01  0.08 
+0.02  0.02 
+0.02  0.05 
+0.02  0.08 


+  NEWAXES  CARTESIAN  . 

+0.02  0.0  0.0  0.0  0.0  0.0. 

+  NODES  . 

+  0.02  0.02  . 

+  0.02  0.08  . 

+0.04  0.02  . 

+  0.04  0.05  . 

+  0.04  0.08  . 

+  . 

+ 

+ 

Although  this  data  file  Is  correct.  It  Is  quite  evident  that  It 

has  no  specific  usefulness  with  respect  to  this  example. 
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SCALE 
EX  EY  EZ 

If  necessary,  this  entry  defines  scale  factors  (EX,  EY,  EZ)  which 
are  applied  to  the  node  coordinates  Included  In  the  following 
NODES  entry(les),  during  the  description  of  the  mesh.  To  switch 
off  this  scaling,  the  SCALE  entry  has  to  be  used  with  the 
parameter?  0,  0,  0  or  1,  1,  1. 


Test  example.  + 

+  NODES  .... 
+0.00  0.02 
+0.00  0.05 
+0.00  0.08 
+  0.01  0.02 
+0.01  0.08 
+0.02  0.02 
+0.02  0.05 
+0.02  0.08 


+  NEWAXES  CARTESIAN  . 

+  0.02  0.0  0.0  0.0  0.0  0.0  . 

+  SCALE  . 

+  2.0  1.0  1.0  . 

+  NODES  . 

+  0.01  0.02  . 

+0.01  0.08  . 

+0.02  0.02  . 

+0.02  0.05  . 

+0.02  0.08  . 

+  . 

+ 

+ 


This  example  Is  a  simple  combination  of  the  entries  NEWAXES  and 
SCALE.  It  allows  the  modification  of  the  thickness  of  the 
electrode  by  simply  modifying  the  scale  factor  (2.00  In  the  above 


•  ELEMENTS 

RAAAAA  MATERA  NGA 

Nil  N12  N13  .  .  NIP 

N21  N22  N23  .  .  N2P 


NK1  NK2  NK3  .  .  NKP 

.  OUMUY  . 

RBBBBB  MATERB  NG8 


Nil  N12  N13 . N1Q 

N21  N22  N23 . N2Q 

NL1  NL2  NL3 . NLQ 


DUMMY 


.  DUMMY  . 

RZZZZZ  MATERZ  NGZ 
Nil  N12  N13  .  N1R 

N21  N22  N23  .  .  N2R 

•  ••••• 
NM1  M42  NM3  .  NMR 

.  DUMMY  . 

.  DUMMY  . 


This  entry  provides  Information  about  the  elements.  It  is 
compulsory,  its  parameters  are  defined  set  by  set,  each  set 
corresponding  to  successive  elements  of  the  same  type.  The 
various  types  of  elements,  their  topologies  and  their 
corresponding  material  and  geometrical  properties  are  described  In 
chapter  4. 

For  each  set  of  parameters,  the  first  line  contains  the  name  of 
the  element  type,  corresponding  to  RAAAAA,  RBBBBB,  ...,  RZZZZZ  In 
the  above  entry  example,  the  name  of  the  corresponding  material 
data  set,  denoted  MATERx  (see  the  entry  MATERIALS  for  more 
details),  and  the  number  of  the  corresponding  geometrical  property 
set,  denoted  NGx  (see  the  entry  GEOMETRY  for  more  details).  Then, 
each  following  line  contains  the  actual  numbers  of  the  element 
nodes,  given  in  a  compulsory  order  which  corresponds  to  the 
element  topology  (one  or  two  lines  per  element,  depending  upon  the 
number  of  nodes  belonging  to  It,  the  first  line  being  terminated 
by  the  character  &  l  f  a  second  line  Is  needed).  Finally,  the  set 
has  to  be  closed  by  a  dummy  line.  Moreover,  the  whole  entry 
parameter  list  has  also  to  be  closed  by  a  dummy  line.  Thus,  the 
two  last  lines  of  this  entry  are  always  two  dummy  lines. 

From  the  computational  point  of  view,  elements  are  assembled  In 
the  sequential  order  of  this  entry  parameter  list.  Thus,  In  most 
of  the  cases,  a  great  deal  of  attention  has  to  be  given  to  this 
order,  which  has  an  Important  effect  on  the  array  sizes  and  the 
CPU  time  requested  by  the  Job. 


Tost  example. 


+ 

+  ELEMENTS  . 

+  R37521  X51  1  . 

+1  3682457  ... 

+  . 

+  R36210  AU4G  2  . 

+6811  13  7  9  10  12 

+  . 

+  . 

+ 

+ 


The  above  numberings  of  the  nodes  can  be  easily  obtained  using  the 
corresponding  element  descriptions  In  chapter  4.  "X51"  and  "AU4G" 
are  parameters  which  refer  to  the  MATERIALS  entry.  "1"  and  "2"  are 
parameters  which  refer  to  the  GEOMETRY  entry.  Their  meaning  Is 
discussed  in  the  two  following  sections. 
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MATERIALS 
MATER 1 
XII  X12  XI 3 
MATER2 
X21  X22  X23 


.  X2P 


MATERN 

XN11  XN12  XN1 3  XN14  XN15  XN16  & 

Ql 

XNM1  XNM2  XNM3  XNM4  XNM5  XNM6  & 
.  DUMMY  LINE  . 


This  entry  provides  information  about  the  material  properties.  It 
Is  compulsory  and  has  to  be  closed  by  a  dummy  line.  Its  parameters 
are  defined  by  set.  Each  parameter  set  corresponds  to  a  given 
material  and  Is  Identified  by  a  name,  selected  by  the  user  and 
containing  less  than  8  characters  (MATER1 ,  MATER2 ,  ....MATERN  in 
the  above  example).  This  name  links  a  set  of  elements  from  the 
previous  ELEMENTS  entry  to  the  given  material.  The  list  of  the 
requested  properties  Is  provided  in  the  element  description,  in 
chapter  4.  Thus,  as  examples,  this  list  contains  Young's  modulus, 
Poisson's  ratio  and  density  in  the  case  of  a  three-dimensional 
elastic  element,  and  it  contains  the  whole  elastic,  piezoelectric 
and  dielectric  tensors  in  the  case  of  a  piezoelectric  element.  If 
a  line  Is  insufficient  to  describe  the  properties  of  a  given 
material,  it  can  be  extended  using  the  character  &.  The  order  in 
which  the  physical  properties  have  to  be  written  Is  compulsory  and 
Is  given  below. 

.  elast Ic  mater lal : 

E  NU  RO 

where  E  is  Young's  modulus,  NU  is  Poisson's  ratio  and  RO  Is  the 

dens  I ty . 

.  plezoelectr Ic  mater lal : 


0.0 

0.0 

RO 

0.0 

0.0 

0.0 

S11 

SI  2 

SI  3 

SI  4 

SI  5 

SI  6 

S21 

S22 

S23 

S24 

S25 

S26 

S31 

S32 

S33 

S34 

S35 

S36 

S41 

S42 

S43 

S44 

S45 

S46 

S51 

S52 

S53 

S54 

S55 

S56 

S61 

S62 

S63 

S64 

S65 

S66 

dll 

dl  2 

dl  3 

dl  4 

dl  5 

d16 

d21 

d22 

d23 

d24 

d25 

d26 

d31 

d32 

d33 

d34 

d35 

d36 

EP1 1 

EP 1  2 

EP1 3 

0.0 

0.0 

0.0 

EP21 

EP22 

EP23 

0.0 

0.0 

0.0 

EP31 

EP32 

EP33 

0.0 

0.0 

0.0 

whore  [S]  Is  the  constant  electric  field  elastic  tensor,  [d]  Is 
the  piezoelectric  tensor  and  [EP]  Is  the  constant  strain 
dielectric  tensor. 

.  fluid  material : 

COMP  .0  R0 

where  COMP  Is  the  modulus  of  compression  and  RO  Is  the  density. 


The  material  property  values  of  classical  or  more  frequently  used 
materials  can  be  written  In  a  special  file,  named  MATER. STD.  If  It 
Is  the  case,  these  materials  need  not  be  described  In  the 
MATERIALS  entry  and  their  names  can  be  used  directly  as  parameters 
In  the  ELEMENTS  entry.  If  one  of  the  names  belonging  to  the 
MATER. STD  file  is  used  In  the  MATERIALS  entry,  the  new  property 
values  which  are  listed  as  parameters  will  automatically  replace 
the  old  values  In  the  file  and  will  be  used  by  the  code  during  the 
run.  In  this  case,  a  warning  will  be  printed  In  the  listing  file. 
Finally,  If  a  new  material  name  is  defined  in  the  MATERIALS  entry, 
this  name  and  the  corresponding  property  values  will  be 
automatically  added  to  the  MATER. STD  file. 

The  MATER. STD  file  can  be  edited  by  the  user  for  any  modification, 
such  as  adding  a  new  material  or  correcting  old  property  values. 
Its  first  line  contains  a  control  Integer  (generally  80,  which  is 
the  number  of  columns),  the  date  of  the  last  edition  and  Its  total 
number  of  lines.  The  following  lines  contain  the  material  property 
values,  listed  as  defined  above  (1  line  for  elastic  or  fluid 
materials,  13  lines  for  piezoelectric  materials).  Each  of  these 
following  lines  starts  with  two  special  fields,  containing 
respectively  the  material  name  and  the  number  of  lines  set  apart 
for  this  material  (as  examples:  STEEL  1  or  CERAMIC  13  ...).  If 
this  file  Is  edited,  the  user  has  to  take  care  of  the  final  total 
number  of  lines  and,  if  needed,  to  correct  its  value  on  the  first 
line  of  the  file.  An  opening  error  can  take  place  during  the 
following  run  of  ATILA  If  this  number  is  wrong. 


15 


Test  example. 


+ 

+  MATERIALS  . 

+  AU4G  . 

+  .714E+11  .344  2780 . 

+  X51  . 

+  0.0  0.0  7350.  0.0  0.0  0.0  . 

+  0.114E-10  -0 . 339E-1 1  -0.410E-11  0.0  0.0  0.0  &  . 

+  -0 . 339E-1 1  0.114E-10  -0.410E-11  0.0  0.0  0.0  &  . 

+  -0.410E-11  -0.410E-11  0.126E-10  0.0  0.0  0.0  &  . 

+  0.0  0.0  0.0  0.510E-10  0.0  0.0  &  . 

+  0.0  0.0  0.0  0.0  0.510E-10  0.0  &  . 

+  0.0  0.0  0.0  0.0  0.0  0.296E-10  8.  . 

+  0.0  0.0  0.0  0.0  0.491E-09  0.0  &  . 

+  0.0  0.0  0.0  0.491E-09  0.0  0.0  &  . 

+  -0 . 955E-1 0  -0 . 955E-1 0  0.208E-09  0.0  0.0  0.0  &  . 

+  0 . 667E-08  0.0  0.0  0.0  0.0  0.0  &  . 

+  0.0  0.667E-08  0.0  0.0  0.0  0.0  &  . 

+  0.0  0.0  0.687E-08  0.0  0.0  0.0  . 


GEOMETRY 

NG1 

XII  X12 
NG2 

X21  X22  .  .  X2P 


NGM 

XM1  XM2  .  .  XMQ 

_  DUMMY  LINE  . 


This  entry  provides  Information  about  the  geometrical  properties. 
It  must  be  closed  by  a  dummy  line.  Its  parameters  are  defined  set 
by  set,  each  parameter  set  corresponding  to  a  given  geometry  and 
being  Identified  by  a  number  selected  by  the  user  (NG1  , 
NG2 .... , NGM  In  the  above  example).  For  a  given  geometry,  the  list 
of  the  requested  properties  Is  provided  In  the  element 
description,  In  chapter  4.  Thus,  as  examples,  this  list  contains 
the  thickness  In  the  case  of  a  plate  element,  and  It  contains  the 
radius  of  curvature  in  the  case  of  a  damping  element.  The  order  In 
which  the  geometrical  properties  have  to  be  written  Is  compulsory. 


Test  example. 


In  this  case,  the  correct  entry  for  geometrical  properties  Is  the 
GEOMETRY  POLARIZATION  entry  which  Is  described  In  the  next 
section.  In  fact,  its  format  Is  Identical  to  the  GEOMETRY  entry 
format,  and  thus  the  example  can  be  considered  as  common  to  the 
two  entr ies. 


GEOMETRY  POLARIZATION  STRING 
NG1 

XII  XI 2 
NG2 

X21  X22  .  .  X2P 


NGM 

XM1  XM2  .  .  XMQ 

_ DUJJ4Y  LINE  . 


This  entry  provides  information  about  the  geometrical  properties 
when  piezoelectric  elements  are  present  In  the  mesh,  and  It  is 
compulsory  in  this  case.  The  GEOMETRY  entry  features  are  all 
relevant  to  this  new  entry  and  must  be  used  for  non-p I ezoe I ectr  I c 
elements.  However,  for  piezoelectric  elements,  specific 
parameters  have  to  be  added  to  define  the  polarization  type 
(cartesian,  cylindrical  or  spherl'cal)  and  or  ientat  ion. The 
polarization  type  Is  specified  by  replacing  STRING  by  CARTESIAN, 
CYLINDRICAL  or  SPHERICAL  in  the  entry  name. 


.  If  STRING-CARTESIAN,  the  polarization  is  homogeneous  In  the 
whole  element  and  the  natural  axes  0X-( ,  OX2  and  0X3  have  to  be 
specified  with  respect  to  the  global  axes  OX,  OY  and  OZ  (0X3  is 
the  polarization  axis  1.1  the  classical  case  of  a  ceramic 
material).  To  do  this,  the  user  has  to  provide  three  Euler  angles, 
named  ALPHA,  BETA  and  GAMMA,  which  transform  OXYZ  into  OX3XiX2  and 
are  defined  as  follows  (see  figure  1): 

-  ALPHA  Is  the  angle  of  the  first  rotation,  around  OZ ,  which 
transforms  the  OXYZ  system  into  a  OX'Y'Z  system  such  that  OX' 
belongs  to  the  OZX3  plane, 

-  BETA  is  the  angle  of  the  second  rotation,  around  OY',  which 
transforms  the  OX'Y'Z  system  Into  a  OX3Y'Z’, 

-  GAMMA  is  the  angle  of  the  third  rotation,  around  0X3,  which 
transforms  the  OX3Y'Z’  system  Into  the  OX1X2X3  system.  Its  value 
Is  dummy  in  the  particular  case  of  a  ceramic  material  (6mm 
symmetry  c  I  ass) . 

The  entry  parameters  are  ALPHA,  BETA,  GAMMA,  listed  in  this 
compulsory  order.  Attention  has  to  be  given  to  the  exact  meaning 
of  this  "cartesian"  polarization.  In  fact.  If  the  polarization  Is 
homogeneous  In  a  three-dimensional  element,  it  is  physically 
homogeneous  In  the  corresponding  part  of  the  material.  On  the 
other  hand,  If  the  polarization  Is  homogeneous  and  Its  direction 
orthogonal  to  the  symmetry  axis  In  the  case  of  an  ax  I  symmetr  I  ca  I 
element,  this  polarization  Is  physically  cylindrical  In  the 
corresponding  part  of  the  material. 


1  8 
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.If  STRING-CYLINDRICAL,  the  polarization  Is  radial  around  a  given 
axis  O'Z*  In  the  whole  element.  The  user  has  first  to  specify  a 
O'X'Y'Z*  system  which  contains  the  O'Z'  axis.  To  do  this,  he 
provides  the  0'  coordinates  (XO,  YO,  ZO)  as  well  as  the  Euler 
angles  (ALPHA,  BETA,  GAMMA)  which  transform  the  OXYZ  system  Into 
the  O'X'Y'Z'  system  and  are  defined  as  follows  (see  figure  2a): 

-  ALPHA  Is  the  angle  of  the  first  rotation,  around  OZ ,  which 

transforms  the  OXYZ  system  into  a  OxyZ  system  such  that  Ox 
belongs  to  the  OZX'  plane, 

-  BETA  is  the  angle  of  the  second  rotation,  around  Oy ,  which 

transforms  the  OxyZ  system  Into  a  OX'yz, 

-  GAMMA  Is  the  angle  of  the  third  rotation,  around  OX',  which 

transforms  the  OX'yz  system  into  the  OX'Y'Z'  system. 

It  Is  self  evident  that  the  choice  of  the  O'X'  (or  O'Y')  axis  is 
arbitrary  and  has  to  be  the  simplest.  With  these  new  reference 

axes,  the  angle  of  the  rotation  around  O'Z'  which  brings  the  O'X' 
axis  to  coincidence  with  the  natural  MX3  axis,  for  each  point  M 
which  is  concerned.  Is  directly  computed  by  the  code  (see  figure 
2b).  Finally,  the  user  must  also  provide  the  angle  DELTA  which 
brings,  for  every  point  M,  the  MZ'  axis  Into  coincidence  with  the 
MX2  natural  axis  (see  figure  2c).  The  value  of  DELTA  Is  dummy  in 
the  case  of  a  ceramic  material  (6mm  symmetry  class).  The  entry 
parameters  are  ALPHA,  BETA,  GAMMA,  XO,  YO,  ZO,  DELTA,  listed  In 
this  compulsory  order.  Attention  has  to  be  given  to  the  exact 
meaning  of  this  “cy I  I ndr I ca I"  polarization.  In  fact,  If  the 
polarization  is  cylindrical  in  a  three-dimensional  element,  It  Is 
physically  cylindrical  In  the  corresponding  part  of  the  materia1. 
On  the  other  hand,  If  the  polarization  Is  cylindrical  around  an 
axis  normal  'to  the  plane  of  an  ax Isymmetr lea  I  element  and 
Intersecting  the  symmetry  axis  at  the  point  0,  this  polarization 
Is  physically  spherical  around  0  In  the  corresponding  part  of  the 
material. 


.  if  STRING-SPHERICAL,  the  polarization  is  radial  around  a  given 
point  O'  In  the  whole  element.  The  user  has  first  to  specify  a 
O'X'Y'Z'  system  centered  on  the  origin  O'.  To  do  this,  he  provides 
the  0'  coordinates  (XO,  YO,  ZO)  and  the  Euler  angles  which 
transform  the  OXYZ  system  Into  the  O'X'Y'Z’  system  (ALPHA,  BETA, 
GAMMA)  and  are  defined  as  for  the  cylindrical  case  (see  figure 
3a).  it  is  self  evident  that  the  choice  of  these  new  axes  is 
arbitrary  and  has  to  be  the  simplest.  With  these  new  reference 
axes,  the  angles  of  the  rotations  which  bring  the  O'X'  axis  into 
coincidence  with  the  natural  MX3  axis,  for  each  point  M  which  is 
concerned,  are  computed  by  the  code  (see  figure  2b).  Finally,  the 
user  has  also  to  provide  the  angle  DELTA  which  brings,  for  every 
point  M,  the  MZ '  axis  Into  coincidence  with  the  natural  MX2  axis 
(see  figure  2c).  The  value  of  DELTA  Is  dummy  in  the  case  of  a 
ceramic  material  (6mm  symmetry  class).  The  entry  parameters  are 
ALPHA,  BETA,  GAMMA,  XO.  YO,  ZO.  DELTA,  listed  In  this  compulsory 
order  . 


Test  example. 


+ 

+  GEOMETRY  POLARIZATION  CARTESIAN 

+  I  . 

♦  0.0  0.0  0.0  0.0  0.0  0.0  0.0  . .  . 

+  2  . 

+  1 . 

+  . 

+ 

♦ 


As  can  be  seen  after  "I",  the  three  Euler  angles  are  equal  to 
zero,  since  the  polarization  direction  and  the  global  OX  axis  are 
Identical.  Moreover,  the  value  of  the  thickness  after  "2"  Is 
dummy,  since  the  model  Is  ax  1 symmetr I ca I . 


REDUCTIONS 
I  N  MDF 


This  entry  defines  the  frontal  reduction  when  master  degrees  of 
freedom  are  requested  (modal  and  harmonic  analysis).  In  this  case. 
It  Is  compulsory. 

.  |  is  a  control  Integer.  If  1-0,  backsubst i tut  Ion  Is  processed  to 
provide  NLO  full  eigenvectors  from  the  eigenmode  number  N.  If  1-1, 
backsubst i tut  ion  Is  processed  to  provide  NLO  full  eigenvectors  the 
e  I  genf  requenc  I  es  of  which  are  above  N  (In  this  case,  N  Is  an 
integer  value  of  a  threshold  frequency,  given  In  Hz).  The  control 
Integer  N  is  mainly  devoted  to  the  elimination  of  the  rigid  body 
modes,  which  are  provided  by  the  eigenmode  computation,  from  the 
backsubst  l  tut  Ion .  In  some  cases,  It  can  be  used  to  restrict  the 
backsubst I  tut  Ion  to  a  specific  given  mode.  Attention  has  to  be 
paid  to  the  fact  that  only  the  twenty  first  eigenvectors  can  be 
backsubst l tuted ,  Including  the  rigid  body  modes. 

.  MDF  Is  the  chosen  number  of  master  degrees  of  freedom  which  are 
automatically  selected  by  the  code.  Attention  has  to  be  paid  to 
the  fact  that  electrical  potential  degrees  of  freedom  are  always 
condensed,  except  the  applied  potential  In  the  case  of  an  harmonic 
analysis  which  Is  a  master  degree  forced  by  the  code.  Moreover, 
manually  selected  degrees  of  freedom  can  be  defined  In  the  set  of 
data  which  specifies  the  boundary  conditions.  Finally,  the  degrees 
of  freedom  corresponding  to  damping  elements  are  all  master 
degrees  and  are  forced  by  the  code.  Then,  the  sum  of  the  total 
master  degree  number  IDF  (automatic  masters  +  applied  potential 
degree  +  degrees  from  damping  elements)  and  the  automatically 
reduced  electrical  degree  number  must  be  lower  than  the  total 
degree  of  freedom  number  IDT. 


Test  example.  + 

+  REDUCTIONS  . 

+  0  1  23  . 

+ 

+ 

Thus,  all  the  available  degrees  of  freedom  are  master  degrees  and, 
due  to  the  fact  that  no  rigid  body  mode  exists  In  this  case,  the 
five  first  full  eigenvectors  will  be  provided. 


FREQUENCIES 

FI  F2  F3  . FN 

This  entry  Is  only  requested  for  an  harmonic  analysis  but  then  Is 
compulsory.  Its  parameters  are  the  values  of  the  frequency  for 
which  a  computation  Is  needed.  The  maximum  number  of  these  values 
Is  20  If  there  Is  no  damping  (in-alr  harmonic  analysis)  and  10  If 
damping  Is  taken  Into  account  ( In-water  harmonic  analysis).  The 
parameter  line  can  be  extended  using  the  character  &. 


Test  example. 

There  Is  no  test  example  for  this  entry  since  this  entry  Is  only 
concerned  with  an  harmonic  analysis. 


RADIATION  MONOPOLAR  or  DIPOLAR 

This  entry  Is  only  requested  for  the  modelling  of  a  radiating 
structure  but  then  Is  compulsory. 


Test  example. 

There  Is  no  test  example  for  this  entry  since  this  entry  Is  only 
concerned  with  an  harmonic  analysis  of  a  radiating  structure. 


PRINTING 

NPR 

This  entry  defines  the  printing  level  NPR,  which  must  equal  to  1, 
2  or  3,  and  corresponds  to  the  amount  of  Information  available  on 
the  result  listing.  If  NPR-1 ,  only  the  computation  results  are 
printed.  If  NPR-2,  the  more  important  data  are  recalled  and 
printed  together  with  the  computation  results.  This  level  Is  the 
normal  (default)  level.  NPR-3  corresponds  to  the  trace  level  and 
is  only  used  to  detect  faults. 


Test  example. 

For  this  example,  the  default  case  Is  the  best  choice. 


•  END 


\  r  Oescr  rp  VI  on  of  a  da  fa  t  IT e . 


This  entry  marks  the  end  of  the  entry  list  and  Is  compulsory. 


Test  example. 


+ 

+  END 
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3.4  THE  LOADING  AND  BOUNDARY  CONDITION  DATA 


The  data  which  are  associated  with  the  loading  and  boundary 
conditions  and  which  are  described  hereafter  have  to  be  written  In 
fixed  formats.  Work  is  In  progress  to  provide  these  conditions 
using  specific  entries  In  the  next  ATILA  version.  In  all  the  data 
lines,,  a  given  degree  of  freedom  Is  marked  by  the  number  D,  with 
the  following  correspondence:  1  for  UX  (or  P),  2  for  UY,  3  for  UZ , 
4  for  NX  (or  V),  5  for  NY  and  6  for  NZ .  Moreover,  lines  which  are 
parallel  to  a  global  coordinate  axis  and  planes  which  are  normal 
to  a  global  coordinate  axis  are  also  marked  by  a  number,  denoted 
P.  Then,  the  correspondence  is:  1  for  a  plane  normal  to  the  OX 
axis,  2  for  a  plane  normal  to  the  OY  axis,  3  for  a  plane  normal  to 
the  OZ  axis,  4  for  a  line  parallel  to  the  OX  axis,  5  for  a  line 
parallel  to  the  OY  axis  and  6  for  a  line  parallel  to  the  OZ  axis. 


3.4.1  DESCRIPTION  OF  THE  LOADING  CONDITION  DATA. 


In  the  case  of  a  static  analysis,  the  lumped  applied  forces  (or 
moments)  have  to  be  described  using  one  of  the  two  following 
procedures.  The  procedure  Is  selected  using  a  control  number, 
named  CONT,  which  has  to  be  written  first.  In  F10.3  format: 

. CONT . 


I 


Then,  If  CONT-7777.0,  the  nodes  where  the  forces  (or  moments)  are 
applied  are  defined  by  their  coordinates  and,  for  each  applied 

force  (or  moment),  the  following  line  has  to  be  written  (3F10.3, 
215,  F10.3): 

. X . Y . Z D..NLO . F . 

X,  Y,  Z  are  the  coordinates  of  the  fulcrum,  D  is  the  application 
direction  (1  to  3  for  an  applied  force,  4  to  6  for  an  applied 

moment,  following  the  above  mentioned  correspondence  rules),  NLO 
is  the  loading  case  number  (from  1  to  the  value  defined  previously 
by  the  NLOAD  entry),  F  Is  the  force  (or  moment)  value. 

If  CONT-8888.0,  the  nodes  where  the  forces  (or  moments)  are 

applied  are  defined  by  their  number  and,  for  each  applied  force 
(or  moment),  the  following  line  has  to  be  Issued  (315,  F10.3): 

.NODE _ D.  .NLO . F . 

NODE  Is  the  application  node  number,  D,  NLO  and  F  have  the  same 
meanings  as  In  the  previous  case. 

In  each  case,  the  set  of  loading  condition  data  lines  has  to  be 
closed  by  a  dummy  line. 
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WARNING 


If  the  model  Is  ax Isymmetr leal ,  the  force 
(or  moment)  actual  value  has  to  be  divided 
by  2*P I . 

In  the  case  of  a  distributed  loading,  the 
lumped  nodal  forces  (or  moments)  have  to  be 
computed  using  the  corresponding  element 
shape  functions. 


Test  example. 

The  test  example  Is  a  modal  analysis  and  does  not  require  any 
loading  condlt Ion. 


3.4.2  DESCRIPTION  OF  THE  BOUNDARY  CONDITION  DATA. 


The  boundary  conditions  are  generally  associated  with 
mechanical  clampings,  symmetry  planes,  electrodes,  pressure 
released  surfaces...  They  are  defined  by  specific  data  lines,  the 
classical  format  of  which  Is  (315,  4X,  A1 ) .  The  available  cases 
are  arranged  In  two  groups.  The  first  group  contains  only  boundary 
conditions  which  refer  to  the  global  coordinate  axes.  The  second 
group  contains  boundary  conditions  which  require  specific  local 
coordinates. 


Boundary  conditions  defined  In  global  axes. 

*  The  degree  of  freedom  D  Is  deleted  at  the  node  N. 

_ N _ D . 

*  The  degree  of  freedom  D  Is  deleted  for  all  the  nodes  which 
belong  to  a  specific  plane  or  line  P  containing  the  node  N. 

. .  .-N _ D _ P . 

This  condition  can  be  restricted  to  the  solid  or  the  fluid  part  of 
a  mesh  In  the  case  of  coupled  problems  by  specifying  S  or  F  In  the 
fourth  field  and  using  one  of  the  two  following  lines: 

.  .  .-N _ D _ P _ S . 

.  .  .-N _ D _ P _ F . 

*  The  degree  of  freedom  D  Is  deleted  for  all  the  nodes  In  the 
mesh . 


D 


t  f  / 
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*  Degrees  of  freedom  of  type  D  are  Identical  for  all  the  nodes 
which  belong  to  a  specific  line  or  plane  P  containing  the  node  N. 


—  -N _ D -P . 


*  The  degree  of  freedom  D  at  node  N  Is  a  master  degree  of  freedom. 


_ N. . ,-D. 


*  Degrees  of  freedom  of  type  D  are  master  degrees  of  freedom  for 
all  the  nodes  which  belong  to  a  specific  line  or  plane  P 
containing  the  node  N. 


.-N. . .-D. 


*  Degrees  of  freedom  of  type  D  are  master  degrees  of  freedom  and 
Identical  for  all  the  nodes  which  belong  to  a  specific  line  or 
plane  P  containing  the  node  N. 


. . .-N — -D. . .-P. 


*  In  an  harmonic  analysis  of  a  piezoelectric  structure,  the  degree 
of  freedom  associated  with  the  applied  potential  has  to  be  a 
master  degree  of  freedom  and  specifically  marked,  since  Its 
numerical  processing  requires  a  special  storage.  This  selection  Is 
realised  using  the  following  data  line: 


•  4  «  ■  N  4  4  4  "“4  • 


where  N  is  a  node  which  belongs  to  the  related  electrode. 
Moreover,  this  marked  degree  of  freedom  Is  the  same  for  the  whole 
electrode.  Thus,  If  this  electrode  Is  a  specific  line  or  plane  P, 
the  data  line  can  be  written: 


*  In  an  harmonic  analysis  of  a  piezoelectric  structure  using  the 
mixed  plane  wave-finite  element  method,  the  degree  of  freedom 
associated  with  the  displacement  D  normal  to  the  Interface  between 
the  two  domains  has  to  be  a  master  degree  of  freedom  and 
specifically  marked,  since  Its  numerical  processing  requires  a 
special  storage.  This  selection  Is  realised  using  the  following 
data  line: 


where  N  Is  the  node  which  corresponds  to  D.  Moreover,  this  marked 
degree  of  freedom  is  the  same  for  the  whole  Interface.  Thus,  If 
this  Interface  Is  a  specific  line  or  plane  P,  the  data  line  can  be 
wr l tten : 

4  4  4  -N  .  4  4  -D  4  4  4  -P _ I . 


•  .  *  A  «  -  -  •  - ' 


29 


fcriapter  3.  'Description  of  a  data  f  I  la. 


*  The  degree  of  freedom  of  type  D  at  node  N  Is  Identical  to  the 
degree  of  freedom  of  the  same  type  D  at  node  M. 


Boundary  conditions  defined  In  local  axes. 


Local  axes  can  be  requested  In  two  different  cases.  Firstly, 
the  displacement  which  Is  constrained  can  be  In  a  direction  which 
Is  not  parallel  to  the  global  axis  (figure  1).  Secondly,  -a 
boundary  condition  has  to  be  applied  on  a  line  which  Is  not 
parallel  to  a  global  axis  or  on  a  plane  which  is  not  perpendicular 
to  a  global  axis  (figure  2).  In  fact,  these  two  cases  are 
generally  merged  (figure  3). 


Figure  1 


Figure  2. 


Figure  3. 
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*  To  define  local  axes  at  a  given  node  N ,  the  user  has  to  Issue 
the  following  lines  (formats  215  and  6F10.3). 


_ N.  . .10. 


.A22 - / 

.  .  . A23 


All,  A12,  A13  are  the  direction  cosines  of  the  first  local  axis 
and  A2 i ,  A22,  A23  those  of  the  second  local  axis.  Then,  degrees  of 
freedom  at  node  N  are  automatically  defined  by  the  code  In  the  new 
local  axes  and  they  can  be  constrained  using  the  same  data  lines 
as  for  the  first  group.  If  the  same  global  axes  have  to  be  defined 
for  all  the  nodes  which  belong  to  the  same  line  or  plane  P 
containing  the  node  N,  the  preceedlng  data  set  has  to  be  modified 
by  simply  substituting  for  the  first  line  the  following  line: 


.  .  .-N.  .  .  10 _ P 


*  To  constrain  nodes  belonging  to  a  plane  which  is  not 
perpendicular  to  a  global  axis,  the  corresponding  line  of  the 
first  group  has  to  be  modified  by  substituting  9001  to  P  and 
adding  Immediately  after  a  second  line  which  contains  the 
direction  cosines  A1  ,  A2  and  A3  of  the  normal  to  this  r.ew  plane 
(format  3F10.3).  Thus,  If  the  degrees  of  freedom  of  type  D  of  all 
the  nodes  which  belong  to  this  plane  have  to  be  deleted,  the  two 
data  II nes  are : 

.  .  .-N - D.9001 . 

. A1 . A2 . A3 . 

To  constrain  nodes  belonging  to  a  line  which  Is  not  parallel  to  a 
global  axis,  the  corresponding  line  of  the  first  group  has  to  be 
modified  by  substituting  9004  to  P  and  adding  Immediately  after  a 
second  line  which  contains  the  direction  cosines  A1 ,  A2  and  A3  of 
this  new  line  (format  3F10.3).  The  preceedlng  example  can  be 
Immedlatly  transposed. 

*  Finally,  the  two  preceedlng  cases  can  be  merged  simply  by  mixing 
the  corresponding  data  lines.  Thus.  If  the  preceedlng  local  axes 
have  to  be  defined  for  all  the  nodes  belonging  to  the  preceedlng 
plane,  the  following  data  lines  have  to  be  Issued: 

. . .-N. . .10.9001  . 


All . A12 . A13 . A21 . A22.  .  .  ./ 

.  .  .  A23 

.  A 1 . A2 . A3 . 


General  comments. 


.  This  data  line  set  has  to  be  closed  by  a  dummy  line. 

.  Constraints  on  different  degrees  of  freedom  can  be  merged  on  the 
same  data  line  by  placing  side  by  side  the  corresponding  Integers 
In  the  definition  of  the  D  Integer.  Thus,  If  the  same  constraint 
has  to  be  applied  to  UX,  UY  and  NZ,  D  has  to  be  equal  to  126. 

.  Null  flux  condition  is  the  natural  condition  of  the  finite 
element  formulation  used  in  ATILA  and  Is  Implicit.  This  Is  the 
case  for  the  fluid  velocity  In  fluid  domains  and  for  the  electric 
field  In  piezoelectric  domains. 


Test  example.  + 

+ . 3 . 

+  .  .  .-1  .  .  .14 _ 5 . 

+  .  .  .-6 _ 4 _ 5 . 

+ . 

+ 

+ 

The  first  line  deletes  the  UZ  component  of  the  displacement  field 
(compulsory,  due  to  the  fact  than  LCPDDC-7).  The  second  line  Is 
the  symmetry  condition  along  the  OY  axis.  The  third  line  Is  the 
open  electrode  condition  on  the  upper  side  of  the  ceramic  ring. 


3.5  THE  MIXED  MOOEL  SPECIFIC  DATA 


When  a  mixed  finite  e l ement -p I ane  wave  model  is  used,  the 
finite  element  part  is  described  simply  following  the  guide  imes 
given  In  the  other  sections.  This  section  is  devoted  to  the 
description  of  the  plane  wave  part  of  the  structure,  which  is 
assumed  to  be  organised  as  shown  below. 


/////////// 

/////////// 
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if  I  NO— 1 ,  the  structure  Is  ax l symmetr I ca I .  If  IND-2,  the  structure 
is  three-dimensional  with  one  symmetry  plane.  If  I  NO-3,  the 
structure  Is  three-dimensional  with  two  symmetry  planes.  The 
following  lines  describe  the  successive  slices  of  part  1,  2  and  3. 
Their  format  Is  (215,  6E10.3).  They  have  to  be  written  under  the 
form : 

■  ITO.  .  ITH . TH . 0R1 . IR1 . 0R2 _ / 

.  . .  IR2 . PHI 

where : 

ITO  -  Ml,  where  N  and  M  are  two  Integers 
def i ned  hereafter  , 

ITH  -  0  or  180  and  ai lows  to  reverse  the 
ceramic  polarization  direction,  its  value  is 
dummy  for  other  materials. 


TH  is  the  thickness  of  the  slice. 


OR  1 
face , 

1  s 

t  he 

outside 

radius 

of 

the  slice 

first 

1  R 1 

face  . 

1  s 
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inside 

radius 

of 

the  slice 

first 

OR2  i 
face , 
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outside 

radius 

of 

the  slice 

second 

IR2 

face  . 

1  3 
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inside 

radius 

of 

the  slice 

second 

PHI 

sign; 

l  s 

the  app i i ed 
i  f  any 

potential 

(beware  of  the 

N  def ' nes  the  type  of  the  si  ice 

N  -  1  for  an  elastic  cylinder  (part  1  or  3), 

-  2  for  a  p i ezoe i ect r t c  ceramic  cylinder 
(part  n  , 

-  3  for  an  elastic  cone  (part  1  or  3), 

«  -1  -2  or  -3  if  this  slice  is  the  last 

s  ce  n  part  3  , 

-  4  for  an  e'astic  cylinder  (part  2), 

-  5  for  a  conical  cy'mder  (part  2  ) 

-  -4  or  -5  if  this  s i 1 ce  is  the  last  slice  in 
part  2  . 
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Chapter  3.  Description  of  a  data  file 


If  IP3  Is  equal  to  zero,  the  normal  displacement  to  P3  Is  free, 
while  If  IP3  Is  equal  to  1  this  normal  displacement  Is  clamped. 


WARNING:  With  this  version,  only  39  successive  slices  can  be 
accomodated  by  the  code. 


*•  Cylindrical  elements  can  be  replaced  by  elements  with 
different  cross  section  shapes  (squares,  rectangles  ...),  If  the 
plane  wave  assumption  Is  still  correct  in  these  cases.  Then,  the 
Inside  and  outside  radius  of  the  corresponding  cylindrical  slices, 
as  defined  in  page  34,  have  to  be  computed  by  the  user  to  ensure 
that  the  cross  section  area  values  are  the  actual  values. 

**•  Material  property  values  are  provided  by  two  DATA 
Instructions  in  a  subroutine  named  TAC.  The  user  can  modify  them 
by  editing  the  corresponding  FORTRAN  file.  For  elastic  materials, 
values  are  listed  In  the  DATA/SM/  Instruction,  with  two  values  for 
each  material  (density  and  Young's  modulus),  successive  materials 
being  put  in  the  M  ascending  order.  For  piezoelectric  materials, 
values  are  listed  in  the  DATA/SP/  instruction,  with  eleven  values 
for  each  material  (density,  constant  electric  field  elastic 
constants  S33,  S13,  S^,  S 1 2  •  S44,  piezoelectric  constants  ^33, 
d31  •  d15  and  constant  strain  dielectric  constants  EP33  AND  EP^). 
Successive  materials  are  also  put  In  the  M  ascending  order. 


3.6  THE  GRAPHIC  DISPLAY  SPECIFIC  DATA 


With  the  ATILA  code,  graphic  displays  can  be  obtained  using 
different  devices,  mainly  plotters  and  workstations.  At  the 
present  time,  results  have  been  obtained  with  BENSON  plotting 
tables,  VERSATEC  Xerox  tables,  Tektronix  4014,  4114,  4115,  as  well 
as  MASSCOMP  and  MICROVAX  computers.  When  plotters  have  to  be  used, 
the  ATILA  run  generates  a  file  which  contains  the  required 
information  for  the  plotter  and  this  plotter  is  generally 
activated  In  a  following  Job.  In  this  case,  graphic  display 
specific  data  have  to  be  added  to  the  data  file  and  are  described 
hereafter.  When  a  workstation  is  used,  the  graphic  part  of  the  Job 
Is  interactive  and  no  data  which  are  related  to  the  graphic 
display  have  to  be  put  in  the  data  file.  Nevertheless,  during  the 
run  of  the  Interactive  Job,  questions  are  submitted  to  the  user 
which  correspond  strictly  to  the  following  Information. 

Graphic  displays  of  the  structure  In  Its  rest  position  or  in 
deformed  shapes  can  be  easily  obtained  using  the  following  set  of 
data.  For  each  graphic  display,  two  data  lines  are  requested  with 
the  respective  formats  (4F10.3,  315)  and  (A40,  F10.3).  These  lines 
are: 

. ALPHA . BETA . XLTH . BRTH. . .NP. .NE1 .  .NE2 

. TITLE . SCALE . 


•  ALPHA  and  BETA  are  angles  (degrees)  which  define  the  projection 
of  the  structure  on  the  drawing  plane.  The  structure  being 
Initially  observed  by  the  user  In  the  OY  direction,  out  of  this 
user,  the  OX  axis  Is  rotated  of  the  ALPHA  angle  around  the  OZ 
axis.  Then,  the  Initial  OZ  axis  Is  rotated  of  the  BETA  angle 
around  the  original  OX  axis.  Finally,  the  structure  Is  projected 
In  the  Initial  OXZ  plane.  ALPHA  and  BETA  are  positive  If  the  axis 
are  moving  back  of  the  Initial  OXZ  plane. 

•  XLTH  and  BRTH  (In  meters)  define  the  space  occupied  by  the 
graphic  display.  If  BRTH  equals  zero  and  XLTH  Is  negative,  -XLTH 
is  the  scale  of  the  display  with  respect  to  the  actual  size  of  the 
structure. 

•  NP  defines  the  type  of  display,  with  the  following 
correspondences : 

NP  -  0  .  Exploded  view  of  the  mesh. 


NP  - 

2  .  Exploded  view  of 

the 

mesh , 

with  the 

node 

number  1 ng . 

NP  — 

•2  .  Exploded  view  of 

t  he 

mesh , 

with  the 

element  numbering. 


<v. 


NP  -  3  .  Dashed  view  of  the  mesh,  with  the 
degrees  of  freedom  numbering. 

NP  -  7  .  Simple  display  of  the  mesh  (full 
I Ines) . 

NP  -12  .  Display  of  the  mesh  with  full  outside 
lines  and  dashed  Inside  lines. 

NP  -31  to  49  .  Graphic  display  of  the  deformed 
shape  of  the  structure.  In  the  loading  case, 
full  eigenvector  or  excitation  frequency  order. 

The  deformed  shape  is  In  full  line  while  the 
rest  position  Is  In  dashed  line.  If  damping  Is 
taken  into  account  In  the  computation,  two 
graphic  displays  have  to  be  Issued  for  each 
frequency  (real  and  Imaginary  part  of  the 
displacement  field). 

For  NP  -  0,  2,  -2,  3,  7  and  12  the  display  Is  provided  by  a 
preprocessor.  If  NP  -  31  to  49,  the  display  Is  provided  by  a 
postprocessor.  Several  displays  can  be  simultaneously  requested  in 
each  of  these  two  categories. 

•  NE 1  and  NE2  permit  getting  the  graphic  display  of  a  part  of  the 
structure  only,  in  this  case,  NE1  Is  the  number  of  the  first  drawn 
element  and  NE2  the  number  of  the  last.  If  these  integers  are 
equal  to  zero,  the  whole  stucture  Is  displayed.  Care  must  be  given 
to  the  case  of  fluid  domain,  since  In  this  case  the  finite  element 
variable  is  the  pressure  which  cannot  be  represented  as  a 
displacement.  So,  fluid  elements  have  to  be  excluded  from  graphic 
displays . 

•  A  specific  TITLE  can  be  given  to  each  graphic  display  This 
title  cannot  contain  more  than  40  letters. 

•  SCALE  Is  a  multiplying  factor  for  the  displacement  amplitude  it 
allows  adjustment  of  this  amplitude  value  for  sake  of  clarity  in 
the  case  of  an  eigenmode  computation,  for  which  the  sign  of  the 
whole  displacement  field  is  not  absolutely  defined,  it  allows  also 
to  reversal  of  this  sign  for  the  sake  of  comparison  between 
results  of  different  computations  {  in  the  case  of  piece  part 
analysis,  for  example) 
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Chapter  4.  Element  library 


4.1  INTRODUCTION 


01 


chapter  describes 
the  AT  I  LA  code*.  For 


Th  I  s 
I  n 

name .  its  name  l n 
def i n i t ion ,  the  I 


the  finite  elements  which  are  available 
each  element,  it  provides  Its  present 
the  previous  ATILA  version  (In  parenthesis/,  its 
st  of  the  degrees  of  freedom  which  are  attached 


to  its  nodes  (translations,  rotations,  electrical  potential, 
pressure),  its  entry  parameters  (topology,  material  properties, 
geometrical  properties). 

Most  of  the  elements  in  the  ATILA  library  are  Isoparametric. 
Thus,  complex  structures  with  curved  sides  or  faces  can  be  easily 
modelled  Nevertheless,  best  results  are  obtained  when  these 
elements  have  reasonable  shapes.  Wild  distortions  lead  to 
inaccuracy  because  the  assumptions  made  within  the  code  for 
strains  become  unrealistic.  Extreme  distortion  may  even  cause 
program  failure  The  user  should  guard  against  excessive 
distortion  by  adhering  to  the  following  guide  lines  for  2D 
elements  or  for  faces  of  3D  elements 


the  angles  between  adjacent 
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4.2  AT  I  LA  ELEMENT  DIRECTORY 


NATURE 

GEOM 

TYPE 

N 

NAME 

COMMENT 

n 

ELASTIC 

HEXA 

3D 

20 

HEXA20E 

4 

PRIS 

3D 

15 

PRIS15E 

5 

QUAD 

2D&AX 1 

8 

QUAD08E 

6 

QUAD 

PLATE 

8 

PLAT08E 

QUAD 

SHELL 

8 

FACE08E 

FACET 

TR  IA 

2D&AX 1 

6 

TRIA06E 

TR  IA 

PLATE 

6 

PLAT06E 

10 

LINE 

SHELL 

3 

SHEL03E 

AX  1 

11 

AX  1 

5 

TRAN05E 

TRANS 

12 

LINE 

SPRING 

2 

SPR I02E 

13 

PIEZO. 

HEXA 

3D 

20 

HEXA20P 

n 

PRIS 

3D 

15 

PR IS15P 

if 

QUAD 

AX  1 

8 

AX IS08P 

Di 

QUAD 

2D 

8 

QUAD08P 

PL.  ST 

17 

TR  IA 

AX  1 

6 

AX  1 S06P 

18 

TRIA 

2D 

6 

TRIA06P 

PL.  ST 

19 

FLUID 
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3D 

20 

HEXA20F 

20 

PRIS 

3D 

15 

PRIS15F 

21 

QUAD 

AX  1 

8 

AXIS08F 

22 

TRIA 

AX  1 

6 

AX  1 S06F 

23 

SPECIAL 

QUAD 

3D 

16 

QUAD 16 1 

1 

24 

TRIA 

3D 

12 

TRIA12I 

25 

LINE 

AX  1 

6 

AX  1 S06 1 

INTER 

26 

QUAD 

3D 

8 

QUAD08M 

MONOP 

27 

TRIA 

3D 

6 

TR 1 A06M 

MONOP 

28 

LINE 

AX  1 

3 

AX IS03M 

MONOP 

29 

L  INE 

AX  1 

3 

AX IS03D 

DIPOL 

30 

GEOM 

:  geome  t  r I ca l  shape 

HEXA 

hex  aedr  a . 

of  t he  e l ament . 

PRIS 

t  r  i  angu  tar  prism. 

N 

number  of  nodes 

QUAD 

auadr i lateral . 

NAME 

name  of  the  element. 

TRIA 

t  r  i  ang i e . 

as  required  for  the 

INTER 

inter  face  e lament . 

entry 

MONOP 

monopolar  damper 

P 

page  number 

DIPOL 
PL  .ST 
TRANS 

dipolar  damper  . 
p i ane  strain, 
t  r  ans i t i on  el ement 
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Chapter  4.  Element  library. 


HEXA20E  (R37110). 


DESCRIPTION. 

. HEXA20E  Is  a  twenty  node  Isoparametric  hexaedra,  used  for 
homogeneous  Isotropic  elastic  materials. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  u2 
(3  translat Ions) . 

TOPOLOGY. 


.20  .  .  .  .8 


18  19  . 

15  16 

5  ...  .17  ....  6 


13 
.  10 


3  ...  .12  ....  4 

14 

.  11 


1  ....  9  ....  2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

HEXA20E  MATER 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15  N16  N17  N18  N19 
N20 

where  N1 , . . .N20  are  the  actual  node  numbers  corresponding  to  the 
integers  1....20  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  Is 
the  dens  I ty  (kg/m^) . 

."GEOMETRY"  entry  parameters. 

No . 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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Chapter  4.  Element  library 


PRIS15E  (R37210) . 


DESCRIPTION. 


.PRIS15E  Is  a  fifteen  node  Isoparametric  triangular  prism,  used 
for  homogeneous  Isotropic  elastic  material. 


.The  degrees  of  freedom  are,  for  each  node:  ux ,  uy,  u2 
(3  translations). 


TOPOLOGY . 


15  14 

12 


.13  ...  5 


11 

9  8 


0  .  .  .  X 


.  .  7 


ENTRIES. 


."ELEMENTS"  entry  parameters. 

PRJS15E  MATER 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15 

where  N1....N15  are  the  actual  node  numbers  corresponding  to  the 

Integers  1....15  on  the  above  figure. 


."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young’s  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 
the  dens  I ty  (kg/m^) . 


."GEOMETRY"  entry  parameters. 
No . 


COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3 


A‘  yy  «  j 
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QUAD08E  (R36210). 


DESCRIPTION. 

.QUAD08E  Is  an  eight  node  Isoparametric  quadrilateral,  used  for 
homogeneous  isotropic  elastic  material  and  for  plane  stress,  plane 
strain  or  ax  I symmetr I ca I  analysis. 

.The  degrees  of  freedom  are,  for  each  node:  ux ,  uy 
(2  translat Ions) . 

TOPOLOGY . 


Y 


.  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 
the  dens  I ty  (kg/m3) . 

."GEOMETRY"  entry  parameter. 

T 

where  T  is  the  thickness  of  the  element  (Ignored  in  the  case  of 
Diane  strain  and  ax l symmet r I ca I  analyses). 

COMMENTS. 


I 


.  Gene  rail 


y . 


the  aspect  ratio  has  to  be  less  than  3 


» 
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PLAT08E  ( R 38 1 10). 


DESCRIPTION. 


.PLAT08E  Is  an  eight  node  flat  quadrilateral  plate  element,  used 
for  homogeneous  Isotropic  elastic  material.  It  relies  on  the 
classical  Love-K i rchhof f  hypothesis. 

.The  degrees  of  freedom  are,  for  each  node:  uz,  nx,  ny 
(1  translation,  2  rotations). 

TOPOLOGY . 


3  ....  8  ....  4 
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O  X 
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5 


.  2 


ENTRIES 

."ELEMENTS"  entry  parameters. 

PLAT08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

*nere  N1  .  N8  are  the  actual  node  numbers  corresponding  to  the 
ntegers  '  .  8  on  the  above  figure. 

mateb  a^S'  entry  parameters. 

E  NU  RO 

**here  E  s  'Dung  s  modu  us  (Pa).  NU  is  Poisson's  ratio  and  RO  Is 
*  np  de-'s  t  *  *  g  r3 

jEOme^P”'  ent'v  parameter 
T 

«yn»'e  "  s  t  Ke  t  n  c«ness  o f  the  element. 

COLNE  NTS 


Se^e'a  >.  assert  rat'C  has  to  be  less  than  3. 

'"9  png  t^d  na  a^d  atera  dimensions  of  the  whole  plate  have 
•r  he  a*  east  ten  t  mes  the  thickness 

s  e  eme-*  ca"  be  super  ■ mposed  to  the  QUAD08E  element  (plane 
c'-fss  rc^d  *  b"  *  r  gene'ate  a  'at  shell  element  (facet). 


FACE08E  (R44210). 


DESCRIPTION. 

.FACE08E  Is  an  eight  node  quadrilateral  shell  element,  used  for 
homogeneous  Isotropic  elastic  material.  It  relies  on  the  classical 
Love-K I rchhof f  hypothesis  and  has  to  be  flat  (facet). 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  u2,  nx,  ny 
(3  translations,  2  rotations). 

TOPOLOGY. 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

FACE08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1,...N8  are  the  actual  node  numbers  corresponding  to  the 
integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

w^e^e  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  Is 
•  -e  dens  i  t  y  ( tcg/m^ )  . 

ome'-b  * "  'jntry  parameter. 

<•»■«*  "  •>  *  ne  thickness  of  the  element  . 


.«  «  •  .  • -»  asceot  rat  10  has  to  be  less  than  3. 

-  „  -  .3  -a  a^d  aterai  dimensions  of  the  whole  shell  have 

*  ■  •  * •  Ties  tr*e  thickness. 

-  «...  ■  ,  ,a  ent  to  the  superposition  of  the  QUAD08E 

.  »  •  t  i  on  and  PLAT08E  element. 


ft 


I  J  1  i,»VlllJ,JlV'  >  1  1  * 


»»p 


''.hap**'  4  I  O***'--’ 


TRlAOftt  w  T*- 


OESCR IPT ION 

.  TR I A06E  s  a  s  «  node  s-.u-ann . i  -  -s- : 

homogeneous  i  sot  roo 1  c  e  as’  c  »a'«’  *  a~/J  c  a^e  s*  •  ess  a'  * 
strain  or  a t  i s ymme t r  ■ ca  ana  rs  s 

.The  degrees  o’  fr©eoom  a'e  *or  ear''  ''^de  u,  .  uf 
(2  t  r ans i a  t i ons  ) 

TOPOLOGY . 
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.  .  4  . 
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ENTRIES. 


."ELEMENTS"  entry  parameters. 

TRIA06E  MATER  NG 
N1  N2  N3  N4  N5  N€ 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  is 
the  dens  1 ty  (kg/m3) . 

."GEOMETRY"  entry  parameter. 

T 

where  T  Is  the  thickness  of  the  element  (Ignored  In  the  case  of 
plane  strain  and  ax  I symmetr I ca I  analyses). 

COMMENTS . 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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PLAT  OAE 


Chapter  4  E  I  green  t  I  I  br  ar  / 


DECCR IPT ION 

PLAT06E  •>  a  •>  *  'v;'je  'at  *•  an,ju  ar  pate  e'amert  used  'or 
*  om<yjer>eo>j$  sc  • r  or  c  e  as'  <  mater  a  t  'e  ies  on  the  r  assica 

c  ,e  x  '  hr,. * '  h^'-’  hes  i  s 

’  h«  de-j'  ees  o'  ''eedom  a r  e .  'or  each  node  uz  .  nz  ,  ny 

’  '  ans  a  •  or'  /  '  o  t  a  t  ons 

TOPOLOGY 


3 


Y 


0  •  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

PLAT06E  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
Integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young's  modulus  (Pa),  NU  Is  Poissons  ratio  and  RO  Is 
the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T 

where  T  Is  the  thickness  of  the  element. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  longitudinal  and  lateral  dimensions  of  the  whole  plate  have 
to  be  at  least  ten  times  the  thickness. 

.This  element  can  be  superimposed  to  the  TRIA06E  element  (plane 
stress  condition)  to  generate  a  flat  shell  element  (facet). 


TO 


SHEL03E  (R42123)  . 


DESCRIPTION. 

. SHEL03E  is  a  three  node  ax  I symmet r l ca  I  thin  shell  element,  the 
formulation  of  which  takes  account  of  Its  double  curvature,  it 
relies  on  the  classical  Love-K l rchhof f  hypothesis. 

.The  degrees  of  freedom  are,  for  each  node:  ux ,  uy,  nz 
(2  translations.  1  rotation). 

TOPOLOGY . 


Y 


0  - 


X 


1  ....  2  ....  3 


ENTRIES. 

.  "ELEMENTS"entry  parameters. 

.•  SHEL03E  MATER  NG 

N1  N2  N3 

where  N1 ,  N2,  N3  are  the  actual  node  numbers  corresponding  to  the 
Integers  1,  2,  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  Is 
the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T  R 

where  T  is  the  thickness  of  the  element  and  R  its  radius  of 
curvature  (around  the  axis). 


COMMENTS. 


.The  longitudinal  and  radial  dimensions  of  the  whole  shell  have  to 
be  at  least  ten  times  the  thickness. 


vvvvvvvv-^.rv,v^.v.-^.T-.vr.-,^ 

Chapter  4.  Element  library. 


TRAN05E  ( R43000 ) . 


DESCRIPTION. 


.TRAN05E  Is  a  five  node  ax  I symmet r i ca I  element  which  ensures  the 
transition  between  quadrilateral  or  triangular  ax  I symmet r I ca l 
elements  and  the  Drevlous  thin  shell  ax  I symmetr I ca I  element,  it 
allows  an  accurate  description  of  their  Interaction.  It  includes  a 
solid  part  (nodes  1,  2,  and  3)  and  a  shell  part  (nodes  4  and  5). 


.The  degrees  of  freedom  are,  for  nodes  1,  2  and  3:  ux ,  u 

(2  translations),  and  for  nodes  4  and  5:  u* ,  uy,  nz  (; 

translations,  1  rotation). 


TOPOLOGY . 

Y 

0  .  .  .  X 


ENTRIES. 

. "ELEMENTS"entry  parameters. 

TRAN05E  MATER  NG 
N1  N2  N3  N4  NS 

where  N1  ,  N2  ,  N3,  N4,  N5  are  the  actual  node  numbers 

corresponding  to  the  Integers  1,  2,  3,  4,  5  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  Is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  Is 
the  density  (kg/m^). 

."GEOMETRY"  entry  parameter. 

T  R 

where  T  Is  the  thickness  of  the  element  and  R  Its  radius  of 
curvature  (around  the  axis). 


COMMENTS . 


.The  longitudinal  and  radial  dimensions  of  the  whole  shell  have  to 
be  at  least  ten  times  the  thickness. 


SPRI02E  <  R 30 1 00  )  . 


Chapter  4  Element  library 


DESCRIPTirN. 

.SPRI02E  <s  a  two  node  spring  element,  which  is  only  able  to 
transmit  to  Its  two  limiting  nodes  a  restoring  force  due  to  a 
change  of  its  length,  its  classical  use  is  for  the  modelling  of 
prestress  bolts  in  transducer  design. 

.The  degrees  of  freedom  are,  for  each  node  ux ,  uy,  Uj 
( 3  t  rans 1  at  ions ) . 

TOPOLOGY . 


0  .  .  .  X 


1 . 2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

SPRI02E  MATER  NG 
N1  N2 

where  N1  and  N2  are  the  actual  node  numbers  corresponding  to  the 
integers  1  and  2  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  .0  RO 

where  E  is  Young's  modulus  (Pa)  and  RO  Is  the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

S 

where  S  is  the  cross  section  area. 

COMMENTS. 


.in  ax  I  symmetr  I  ca  I  models,  the  actual  area  S  has  to  be  divided  by 
2*p I .  In  three-dimensional  models,  If  the  rod  belongs  to  N 
symmetry  planes.  Its  actual  area  S  has  to  be  divided  by  N+1 . 

.  TJt>  I  s  element  Is  nothing  more  than  a  classical  spring.  Attention 
has  to  be  paid  to  spurious  modes  which  can  be  due  to  the  simple 


spr I ng  hypothes I s . 


Chapter  4  Eleven)  *  »ry 


HEXA20P  R3’f>' 


DESCRIPTION. 


HEXA20P  is  a  twenty  node  isoparametric  heraedra.  designed  tor  any 
piezoeiectr  ic  mater  iai 


.The  degrees  of  freedom  are,  for  each  node  u„ ,  uy ,  uz ,  v 
(3  translations.  I  electrical  potential) 


TOPOLOGY . 
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ENTRIES. 


."ELEMENTS"  entry  parameters. 

HEXA20P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15  N16  N17  N18  N19 
N20 

where  N1....N20  are  the  actual  node  numbers  corresponding  to  the 
Integers  1....20  on  the  above  figure. 


."MATERIALS"  entry  parameters. 

The  whole  elast Ic, piezoeiectr Ic  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 


."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  Issued  under  the  special  form 
•GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


Chapter  4  t  l 


nt  library 


PR  I  LISP  M’S',' 


DESCRIPT  I QH 

PR  I S 1 5P  s  a  f  1  * i een  'x>de  sopa ' ame  t r  i c  t '  i angu i a  r  pr  ism, 

les  '  gned  for  any  p>e7oe>ect*  °c  tsai#'  'ai 

The  degrees  of  f'eedom  are  for  each  node  Uj  ,  Uy,  u2 ,  v 

3  translations,  i  eiectricai  potential) 

TOPOLOGY .  6 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

PR IS15P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15 

where  N1....N15  are  the  actual  node  numbers  corresponding  to  the 

integers  1....15  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elast ic.plezoelectr lc  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  Issued  under  the  special  form 
"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


15 


AX IS08P  ( R37521  )  . 


DESCRIPT  ION. 

.AXIS08P  is  an  eight  node  isoparametric  quadrilateral,  designed 
for  any  p i ezoe I ec t r  i  c  material  and  used  only  for  ax l symmet r l ca  I 
ana  lysis 

.The  degrees  of  freedom  are,  for  each  node:  ux ,  Uy,  v 
(2  translations,  1  electrical  potential). 

TOPOLOGY . 
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.  X 


3  ....  8  ....  4 


6 


7 


1 


5 


2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS08P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
Integers  1....8  on  the  above  ftgure. 

.“MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  Issued  under  the  special  form 
"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


QUAD08P  (R37501). 


DESCRIPTION. 

.0UAD08P  Is  an  eight  node  Isoparametric  quadrilateral,  designed 
for  any  piezoelectric  material  and  used  only  for  plane  strain 
ana  lysis. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  v 
(2  translations,  1  electrical  potential). 

TOPOLOGY . 


3  ....  8  ....  4 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
Integers  1,...8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  Issued  under  the  special  form 
"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


Chapter  4.  Element  llbrar 


AXIS06P  (R37522)  . 


DESCRIPTION. 

•AXIS06P  Is  a  six  node  Isoparametric  triangle,  designed  for  any 
piezoelectric  material  and  used  only  for  ax Isymmetr lea  I  analysis. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  v 
(2  translations,  1  electrical  potential).  y 

TOPOLOGY . 
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.  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS06P  MATER  NG 
N1  N2  N3  N4  NS  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
Integers  1 , . . .6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY*  has  to  be  Issued  under  the  special  form 
"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


TRIA06P  ( R37502  '  . 


DESCRIPTION. 

.TRIA06P  Is  a  six  node  isoparametric  triangle,  designed  for  any 
piezoelectric  material  and  used  only  for  plane  strain  analysis 

.The  degrees  of  freedom  are,  for  each  node  ut ,  uy .  v 
(2  translations,  i  electrical  potential). 

TOPOLOGY . 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIA06P  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
Integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 
detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS . 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  Issued  under  the  special  form 
"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


HEXA20F  (R37810). 


DESCRIPTION. 

. HEXA20F  is  a  twenty  node  isoparametric  hexaedra,  used  for  fluid 
homogeneous  med  i  a  . 

.The  degree  of  freedom  is.  for  each  node:  p  (the  pressure). 

TOPOLOGY . 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

HEXA20F  MATER 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15  N16  N17  N18  N19 
N20 

where  N1....N20  are  the  actual  node  numDers  cor r espona ' nj  to  f>e 
integers  1 . . . . 20  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  iPa'  ana  RO  s  the 
dens  l  ty  ( kg/m^  )  . 

."GEOMETRY"  entry  parameters 
No. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3 
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PRIS15F  .R37820). 


DESCRIPT  ION. 

PRIS15F  s  a  '  'tee"-  node  sopar  amet  r  i  c  triangular  prism,  used 
*cr  *  -  a  •'cnogeisous  med  a 

.’’'e  degrae  o'  'reedom  s.  'or  each  node  p  ithe  pressure). 
TOPOLOGY 
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ENTR  I  ES 

:  .  :  U:  N '  •  »  ;  * 

PRS15F  WATER 

S'  N2  N3  S4  N5  N6  S'  N8  NS  NIC  S’’  N 1 »  SI  3  SI  4  N15 

• '  e  e  S'  s  *■  a  e  a.  *  ,.a  ••  de  •'..Te-  s  „  •  esrc-d  n$  to  the 

"  '  e  ge  ■  s  n  ;  ■  *  *-e  a  t  "  >  e  ”  j..  ■  e 

ma  ' :  "  ^  **-•'•*  ~  a  ■  a  ’•'e  '  e  ■  s 

COMP  0  RO 

»*'»'  e  .  ,'Wr  s  •  -e  j..  „s  *  ,  ,'T'  fs?  o"  r  a  ar  d  RO  is  the 
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COMMENTS  . 

.General  iv,  the  aspect  rat'O  has  to  be  less  than  3. 
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Chapter  4.  Element  llbrar 


AXIS08F  (R37531 ) . 


DESCRIPTION. 

.AXIS08F  Is  an  eight  node  Isoparametric  quadrilateral,  used  for 
fluid  homogeneous  media  and  for  ax Isymmetr leal  analysis. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY . 


3 


•  8  • 


4 


ENTRIES. 

."ELEMENTS”  entry  parameters. 

AXIS08F  MATER 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
Integers  1 , . . .8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  R0 

where  COMP  Is  the  modulus  of  compression  (Pa)  and  R0  Is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameters. 

No., 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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Chapter  4.  Element  llbrar 


AXIS06F  (R37532) . 


DESCRIPTION. 

.AXIS06F  Is  a  six  node  Isoparametric  triangle  used  for  fluid 
homogeneous  media  and  for  ax Isymmetr leal  analysis. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


3 


Y 


6 
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0 


.  X 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS06F  MATER 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  Is  the  modulus  of  compression  (Pa)  and  RO  Is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameters. 

No,. 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 


QUAD16I  (R37910). 


DESCRIPTION. 

.QUAD16I  is  a  sixteen  node  isoparametr Ic  quadr I latera  I  element 
which  has  to  be  used  to  ensure  the  matching  between  solid  and 
fluid  meshes  along  t'nelr  Interface  In  the  case  of  three- 
dimensional  analyses.  It  Includes  8  solid  nodes  and  8  fluid  nodes 
which  have  the  same  coordinates  as  the  solid  nodes.  It  has  no 
thickness. 

.The  degrees  of  freedom  are,  for  each  solid  node:  ux,  uy,  uz 
(3  translations),  and  for  each  fluid  node:  p  (the  pressure)  . 

TOPOLOGY. 

SOLID  MESH 

7  ...  .16  ....  8 
.  3  ...  .12  ....  4 


14  .  15  . 

2  Y  .10  .11 


..  5  .  .  .  .13  ....  6  . 

0  .  .  .  X  1  ....  9  ....  2 

FLUID  MESH 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD16I 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12  N13  N14  N15  N16 

where  N1....N16  are  the  actual  node  numbers  corresponding  to  the 

integers  1....16  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

No. 

."GEOMETRY"  entry  parameters. 

No. 

COMMENTS. 

.The  nodes  1,  2,  3,  4,  9,  10,  11,  12  are  fluid  nodes.  The  nodes  5, 
6,  7,  8,  13,  14,  15,  16  are  solid  nodes. 

.Attention  must  be  paid  to  the  orientation  of  the  OZ  axis,  from 
the  fluid  domain  to  the  solid  domain. 


TRIA12I  (R37920) . 


DESCRIPTION. 

. TR IA12I  is  a  twelve  node  Isoparametric  triangular  element  which 
has  to  be  used  to  ensure  the  matching  between  solid  and  fluid 
meshes  along  their  Interface  In  the  case  of  three-dimensional 
analyses.  It  Includes  6  solid  nodes  and  6  fluid  nodes  which  have 
the  same  coordinates  as  the  solid  nodes.  It  has  no  thickness. 

.The  degrees  of  freedom  are,  for  each  solid  node:  ux,  uy,  uz 
(3  translations),  and  for  each  fluid  node.-  p  (the  pressure)  . 

TOPOLOGY. 

SOLID  MESH 

6 

.  .3 


12 

Z  Y  .9  .8 

4.  .10 . 5 . 

0  .  .  .  X  1 . 7 . 2 

FLUID  MESH 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIA12I 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  Nil  N12 

where  N1....N12  are  the  actual  node  numbers  corresponding  to  the 
Integers  1 ....  1 2  on  the  above  figure. 

."MATERIALS”  entry  parameters. 

No. 


."GEOMETRY"  entry  parameters. 
No. 


COMMENTS. 


.The  nodes  1,  2,  3,  7,  8,  9  are  fluid  nodes.  The  nodes  4,  5,  6, 
10,  11,  12  are  solid  nodes. 

.Attention  must  be  paid  to  the  orientation  of  the  OZ  axis,  from 
the  fluid  domain  to  the  solid  domain. 
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AXIS06I  (R37541 ) . 


1  >.■ 
o 


DESCRIPTION. 

.AXIS06I  Is  a  six  node  Isoparametric  element  which  has  to  be  used 
to  ensure  the  matching  between  solid  and  fluid  meshes  along  their 
Interface  In  the  case  of  ax isymmetr leal  analyses.  It  Includes  3 
solid  nodes  and  3  fluid  nodes  which  have  the  same  coordinates  as 
the  solid  nodes.  It  has  no  thickness. 

.The  degrees  of  freedom  are,  for  each  solid  node:  ux ,  uy 
(2  translations),  and  for  each  fluid  node:  p  (the  pressure)  . 

TOPOLOGY. 


Y 

SOLID  MESH 

3  ....  6  ....  4 
0...X  1  ....  5  ....  2 

FLUID  MESH 

i 

ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS06I 

N1  N2  N3  N4  NS  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
Integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

No. 

."GEOMETRY"  entry  parameters. 

No. 

COMMENTS. 

.The  nodes  1,  2,  5  are  fluid  nodes.  The  nodes  2,  3,  6  are  solid 
nodes . 

.Attention  must  be  paid  to  the  orientation  of  the  OY  axis,  from 
the  fluid  domain  to  the  solid  domain. 
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Chapter  4.  Element  1 1 bran 


GUAD08M  (R 37930) 


DESCRIPTION. 

.QUAD08M  Is  an  eight  node  Isoparametric  quadrilateral  damping 
element  used  to  take  account  of  a  radiation  condition.  It  Is 
restricted  to  the  monopolar  (spherical)  damping  condition.  It  has 
to  be  spherical  and  attached  to  the  outside  surface  of  a  three- 
dimensional  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 
TOPOLOGY. 


0  .  .  .  X 


1  .  .  .  5  ....  2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08M  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
Integers  1,...8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  Is  the  modulus  of  compression  (Pa)  and  RO  is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  is  the  radius  of  curvature  of  the  element. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.With  respect  to  the  orientation  of  the  whole  mesh,  Its  normal  has 
to  point  outwards. 


TRIA06M  (R37940). 


DESCRIPTION. 

.TRIA06M  Is  a  six  node  Isoparametric  triangular  damping  element 
used  to  take  account  of  a  radiation  condition.  It  Is  restricted  to 
the  monopolar  (spherical)  damping  condition.  It  has  to  be 
spherical  and  attached  to  the  outside  surface  of  a  three- 
dimensional  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 
TOPOLOGY. 
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6  5 

Y 

O...X  1  ....  4  ....  2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIA06M  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS  entry  parameters. 

COMP  .0  RO 

where  COMP  Is  the  modulus  of  compression  (Pa)  and  RO  Is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  Is  the  radius  of  curvature  of  the  element. 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.With  respect  to  the  orientation  of  the  whole  mesh,  Its  normal  has 
to  point  outwards. 
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AXIS03M  (R37551 ) . 


DESCRIPTION. 

.AXIS03M  Is  a  three  node  Isoparametric  damping  element  used  to 
take  account  of  a  radiation  condition.  It  Is  restricted  to  the 
monopolar  (spherical)  damping  condition.  It  has  to  be  circular  and 
attached  to  the  outside  surface  of  an  ax Isymmetr lea  I  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY . 


Y 


0  .  . 


.  X 


1  ....  3  ....  2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS03M  MATER  NG 
N1  N2  N3 

where  N1 ,  N2  and  N3  are  the  actual  node  numbers  corresponding  to 
the  Integers  1,  2  and  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COM P  Is  the  modulus  of  compression  (Pa)  and  RO  Is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  Is  the  radius  of  curvature  of  the  element. 

COMMENTS. 


.To  ensure  the  correct  orientation  of  this  element,  Its  local  OX 
axis  must  be  oriented  In  the  positive  direction  of  the  global  OX 
axis.  Thus,  the  global  x  coordinate  of  node  2  must  be  higher  than 
the  global  x  coordinate  of  node  1  on  the  above  figure. 
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AXIS03D  (T37551 )  . 


DESCRIPTION. 

.AXIS03D  Is  a  three  node  Isoparametric  damping  element  used  to 
take  account  of  a  radiation  condition.  It  includes  the  dipolar 
damping  condition.  It  has  to  be  circular  and  attached  to  the 
outside  surface  of  an  ax Isymmetr lea  I  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


Y 


0 


.  X 


1  ....  3  ....  2 


ENTRIES. 

. "ELEMENTS"entry  parameters. 

AXIS03D  MATER  NG 
N1  N2  N3 

where  N1 ,  N2  and  N3  are  the  actual  node  numbers  corresponding  to 
the  Integers  1,  2  and  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  Is  the  modulus  of  compression  (Pa)  and  RO  is  the 
density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  Is  the  radius  of  curvature  of  the  element. 

COMMENTS. 


.To  ensure  the  correct  orientation  of  this  element.  Its  local  OX 
axis  must  be  oriented  In  the  positive  direction  of  the  global  OX 
axis.  Thus,  the  global  x  coordinate  of  node  2  must  be  higher  than 
the  global  x  coordinate  of  node  1  on  the  above  figure. 
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5.1  INTRODUCTION 


The  general  organisation  of  an  ATILA  Job  Is  described  In 
section  1.3.  This  chapter  provfdes  the  corresponding  organisation 
of  an  ATILA  run  In  the  special  case  of  a  VAX  or  MICROVAX  computer 
using  the  VMS  operating  system.  The  same  organisation  Is  also 
correct  for  other  computer  types  but  with  different  procedures. 
Moreover,  the  part  of 'this  chapter  which  is  devoted  to  the  graphic 
displays  Implies  that  the  GKS  library  Is  available  and  can  be 
linked  during  the  corresponding  part  of  the  run,  the  user  being 
then  working  Interactively  with  the  help  of  a  graphic  terminal 
(Tektronix  4014,  VT240  ...). 


Chapter  S.  Description  of  an  A~ I  LA  run. 


5.2  THE  MAIN  RUN 

4  '  * 


To  ran  an  ATILA  job,  the  user  must  have  written  a  data  file  as 
described  In  chapter  3.  This  data  file  must  be  named  JOBNAME.ATI, 

ATI  being  the  job  type.  Then  the  user  must  run  the  automatic 
program  generator,  named  PGEN,  using  the  entry: 

RUN  PGEN 

During  this  run,  several  questions  are  asked  by  the  code  which 
allow  the  specification  of  the  job  name  (JOBNAME)  ana  the  Job  type 
(ATI),  the  verification  of  the  analysis  type  which  has  been 
Identified  by  the  code,  the  requirement  of  special  data  storages 
for  a  Job  restart  (only  in  the  case  of  an  harmonic  analysis)  cr  a 
postprocessing  (graphic  displays  ...).  PGEN  provides  three  files 
named  JOBNAME. FOR,  JOBNAME.COM  and  JOBNAME. LST.  JOBNAME. FOR 
contains  the  ATILA  main  program.  Generally,  It  can  be  used  without 
any  modification.  Nevertheless,  It  has  to  be  edited  In  few  special 
cases:  the  modal  analysis  of  an  elastic  structure  flooded  In  an 
incompressible  f'uid  (MOD  5),  the  radiated  pressure  computation  • 

for  a  given  displacement  field  (HAR  3)  and  the  use,  by  advanced 
users,  of  new  algorithms  which  are  currently  under  test  (MOD  1, 

MOD  4,  ...).  In  these  cases,  footnotes  are  available  at  the  end  of 
the  corresponding  sections  in  chapter  2  which  provide  all  the 
necessary  information.  JOBNAME.COM  is  used  to  compile  and  link  the 
program.  It  need  not  be  edited.  Finally,  JOBNAME. LST  contains  a 
copy  of  the  data  file,  as  read  during  the  PGEN  run,  as  well  as 
information  about  the  degree  of  freedom  numbering.  It  allows  a 
very  efficient  check  of  the  data  file  and  mainly  of  the  boundary 
conditions.  Generally,  the  user  must  edit  this  list  and  perform 
the  corresponding  check. 

After  the  PGEN  run,  the  user  has  to  compile  the  main  progam  and 
link  It  to  the  requested  libraries.  This  step  Is  simply  realised 
using  the  entry: 

•JOBNAME 

which  provides  the  executable  (J0BNAME.EXE).  Finally,  the  user  has 
to  run  the  program  using  the  entry: 

RUN  JOBNAME 

During  this  run,  the  code  asks  for  the  Job  name  (JOBNAME)  and  the 
Job  type  (ATI)  and  Information  appears  on  the  user's  screen  which 
describes  the  current  status  (number  of  the  element  the  code  Is 
currently  assembling,  name  of  the  solver  the  code  Is  currently 
using  ...)  and  provides  the  timing.  Moreover,  the  code  writes  down 
several  files  named  JOBNAME. LST,  JOBNAME. WRK,  JOBNAME. SY2  and 


Chapter  S.  Description  of  an  AULA  run. 


J0BNAME.SY4.  JOB NAME .L ST  contains  a  partial  copy  of  the  data  file, 
£  depending  upon  the  printing  level  selected  by  the  user  In  this 

S'  data  file,  and  all  the  results  (elgenfrequencies,  displacement  or 

pressure  fields  associated  with  successive  elgenmodes  or 
excitation  frequencies,  electrical  Impedances  ...  ).  The 

JOBNAME . WRK  and  J0BNAME.SY2  files  are  working  arrays  the  code 
always  needs.  They  contain  useful  Information  at  the  end  of  the 
run  only  If  a  restart  has  to  be  performed.  The  JOBNAME. SY4  file 
contains  several  arrays  the  code  needs  for  postprocessing. 
Particularly,  It  can  provide  the  PLO  array  which  contains  the 
successive  displacement  and/or  pressure  fields  as  well  as  the 
CPDDC  array  which  contains  the  node  coordinates  and  the  degree  of 
freedom  numbering.  This  file  can  be  easily  read  by  advanced  users, 
by  editing  the  JOBNAME. FOR  file,  deleting  all  the  subroutine  calls 
and  replacing  the  WRITE  entries  devoted  to  the  SY4  file  by  the 
corresponding  READ  entries. 


5.3  THE  RESTART 


To  restart  an  ATILA  job  Is  often  necessary  when  the  user  is 
performing  an  harmonic  analysis,  since  more  than  10  or  20 
different  frequencies  are  often  required  In  the  frequency  range  of 
Interest.  The  restart  allows  the  computation  of  the  frequency 
response  using  previously  assembled  matrices  and  often  saves  a 
large  amount  of  CPU  time.  Obviously,  It  can  only  be  performed  If  a 
full  run  has  been  realised  before  and  If  the  required  arrays  have 
been  stored.  Then,  the  process  Is  exactly  the  same  as  for  a  main 
run,  except  that  the  user  has  to  select  the  restart  during  the  run 
of  PGEN,  when  the  corresponding  question  appears  on  the  screen. 


5.4  PRE  AND  POSTPROCESSING 


In  this  version  of  the  ATILA  code,  the  available  pre  and 
postprocessors  are  mainly  devoted  to  the  graphic  display  using  the 
GKS  package,  although  many  specific  postprocessors  (farfleld 
computation,  transmitting  voltage  response  and  directivity 
computations  and  displays,  stress  computations  and  displays, 
isovalue  displays  ...)  are  currently  under  test  but  are  not 
documented  for  the  present  time.  At  the  preprocessing  level,  the 
code  can  provide  graphic  displays  of  the  mesh  with  different  types 
of  numbering  (nodes,  elements,  degrees  of  freedom).  At  the 
postprocessing  level.  It  can  provide  the  same  graphic  displays  as 
well  as  distorted  shapes  of  the  structure  which  describe  the 
resulting  displacement  fields.  Detailed  Information  about  these 
graphic  displays  are  given  in  section  3.6.  To  obtain  a  graphic 
display,  the  user  must  make  the  entry: 

RUN  PGRAPH 

/  *  Then,  PGRAPH  reads  the  JOBNAME .AT I  file  and  writes  down  three 

■v-  files  named  D JOBNAME . LST ,  D JOBNAME .FOR  and  DJOBNAME.COM.  Moreover, 

during  the  run,  questions  appear  on  the  screen  which  are  very 
similar  to  the  question  which  appears  during  a  PGEN  run.  The 
D JOBNAME .LST  file  Is  identical  to  the  JOBNAME. LST  file  which  Is 
obtained  after  a  PGEN  run  and  allows  a  detailed  check  of  the  data 
file  if  PGRAPH  Is  run  at  the  preprocessing  level.  The  DJOBNAME.FOR 
file  contains  the  main  program  which  Is  needed  to  obtain  the 
graphic  displays.  It  need  not  be  edited.  Finally,  the  DJOBNAME.COM 
file  Is  used  to  compile  and  link  the  program*.  Generally,  It  need 
not  be  edited. 

After  the  PGRAPH  run,  the  user  must  compile  the  main  program 
and  link  It  to  the  requested  libraries.  This  step  Is  realised 
simply  using  the  entry: 

•DJOBNAME 

which  provides  the  executable  (DJOBNAME.EXE).  Finally,  the  user 
must  run  the  program  using  the  entry: 

RUN  DJOBNAME 

During  this  run,  the  code  asks  the  Job  name  (JOBNAME),  the  Job 
type  (ATI)  and  requires  the  choice  of  the  pre  or  postprocessing 
%  level.  The  preprocessing  level  can  be  used  as  soon  as  the 

nJv's  JOBNAME. ATI  file  Is  available.  The  postprocessing  level  can  only 

'V\V‘  be  used  If  a  full  ATILA  run  has  previously  been  performed  (main 

run  or  restart)  and  If  the  necessary  storages  have  been  made.  The 
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whole  graphic  display  phase  must  be  performed  using  a  graphic 
terminal  and  the  information  about  the  display  type,  Its 
orientation.  Its  title  ....  are  required  during  the  run  by  Vy 

questions  which  are  written  down  on  the  screen. 


*.  The  link  of  the  main  program  to  the  GKS  library  Is  provided 
by  the  file  LNK1D.COM  which  Is  used  by  DJOBNAME.COM.  Before  the 
first  PGRAPH  run,  the  user  must  check  that  this  link  Is  correct, 
taking  account  of  the  GKS  Installation  on  the  corresponding 
computer . 
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I  *  |  |  *dl  *«|  SI 
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Chapter  S.  Description  of  an  AT  I LA  run. 


S.S  THE  AUTOMATIC  MESH  GENERATION 


The  use  of  the  automatic  mesh  generator  MOSAIQUE  Is  described 
in  the  following  chapter.  It  requires  a  JOBNAME.MOS  file  which  Is 
written  down  by  the  user.  Then  the  user  runs  MOSAIQUE  with  the 
entry : 

RUN  MOSAIQUE 

MOSAIQUE  provides  two  files  named  JOBNAME.LST  and  JOBNAME.ATI . 
JOBNAME . LST  contains  Information  which  are  only  needed  If  a 
problem  appears  during  the  generation.  The  JOBNAME.ATI  file  Is  the 
data  file  for  the  ATILA  main  run.  Depending  upon  the  contents  of 
the  JOBNAME.MOS  file.  It  can  be  used  immediately  or  can  be  edited 
by  the  user  to  add  the  last  entries  and  data,  generally  the 
loading  or  boundary  conditions. 
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5.6  SUMMARY 


Automatic  mesh 


Main  run. 


JOBNAME.MOS 
RUN  MOSAIQUE 


JOBNAME.LST 
JOBNAME . AT  I 


JOBNAME . AT  I 

RUN  PGEN 

JOBNAME.LST 
JOBNAME . FOR 
JOBNAME.COM 

•JOBNAME 

JOBNAME . EXE 

RUN  JOBNAME 


JOBNAME.LST 
JOBNAME. WRK 
JOBNAME. SY2 
JOBNAME. SY4 


Restart 


E 

i 


JOBNAME .AT i 
JOBNAME . WRK 
JOBNAME . SY2 
JOBNAME. SY4 

RUN  PGEN 

JOBNAME. LST 
JOBNAME . FOR 
JOBNAME . COM 

•JOBNAME 

JOBNAME . EXE 

RUN  JOBNAME 


JOBNAME. LST  «- 
JOBNAME. WRK  *- 
JOBNAME. SY2  — 
JOB NAME. SY 4  — 


Postprocessing. 


JOBNAME. AT  I 
JOBNAME . SY4 


RUN  PGRAPH 


DJOBNAME 
D JOBNAME 
OJOBNAME.COM 


PH  — 

.  LST  *-~ 
.FOR  *- 
.COM  —I 


•DJOBNAME  ' 

DJOBNAME . EXE 


RUN  DJOBNAME 


Interactive  graphic  display 


Due  to  their  large  size,  the  files  .LST,  .WRK,  .SY2,  .SY4  are  not 
stored  as  successive  versions  If  they  are  provided  several  times 
with  the  same  job  name  but  are  automatically  deleted  by  the  code. 
If,  for  any  reason,  the  user  wishes  to  save  one  of  the  previous 
version  before  a  new  run  with  the  same  job  name,  It  has  to  rename 
the  corresponding  file. 
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Chapter 


Automat le  mesh  generation. 


6.1  INTRODUCTION 


As  mentioned  In  section  1.3.2  and  In  the  Introduction  of  the 
third  chapter,  partial  regularity  of  the  mesh  can  often  allow  the 
use  of  an  automatic  generation,  which  provides  node  coordinates 
and  element  topologies.  In  the  ATILA  code,  this  work  Is  performed 
by  a  specialized  preprocessor  named  MOSAIQUE.  Then,  the  user  has 
only  to  define  a  gross  splitting,  built  with  super  elements,  and 
to  select  their  automatic  splitting  Into  finite  elements 
(Figure  1).  From  the  numerical  point  of  view,  each  super  element 
is  transformed  Into  Its  natural  axes,  using  the  isoparametric 
transformation.  Then,  It  Is  split  Into  finite  elements,  following 
the  corresponding  data  provided  by  the  user,  and  these  finite 
elements  are  transformed  back  to  the  global  axes  using  the  reverse 
Isoparametric  transformation.  Finally,  redundant  nodes  are 
eliminated  by  searching  nodes  which  belong  to  a  same  given  small 
2D  or  3D  window. 

Entries  for  MOSAIQUE  are  provided  by  a  specific  Input  file 
written  by  the  user  and  named  JOBNAME.MOS.  At  the  end  of  the 
MOSAIQUE  run,  the  output  file  which  contains  the  generated  node 
coo-dlnates  and  topologies  Is  named  JOB NAME .AT I  (see  chapter  5). 
Moreover,  entries  and  data  which  are  required  by  ATILA  but  are  not 
necessary  for  the  mesh  generation  can  be  also  read  by  MOSAIQUE,  If 
present  In  the  JOBNAME.MOS  file,  and  be  copied  In  the  JOBNAME .AT I 
output  file.  Thus,  starting  with  a  simple  Input  file 
(JOBNAME.MOS),  MOSAIQUE  is  able  to  provide  the  whole  ATILA  data 
file  (JOBNAME .AT I ) .  Finally,  It  has  to  be  noted  that,  because  the 
number  of  nodes  which  have  to  be  defined  for  MOSAIQUE  Is  often 
small,  this  procedure  allows  the  user  to  modify  very  easily 
geometrical  parameters  (length,  width  ...)  as  well  as  the  mesh 
steps. 

-Most  of  the  entries  and  data  contained  in  the  MOSAIQUE  data 
file  are  identical  to  the  ATILA  entries  and  data.  However,  the 
NODES  and  ELEMENTS  entries  are  modified,  while  two  new  entries  can 
be  used.  These  entries  and  data  are  described  In  section  6.2. 
Moreover,  the  splitting  of  a  super  element  Into  finite  elements 
can  be  obtained  using  different  frames.  The  frames  which  are 
retained  In  MOSAIQUE  are  described  In  section  6.3. 
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6.2  MOSAIQUE  ENTRIES  AND  DATA 


6.2.1.  THE  HEADER. 

Information  contained  In  section  3.2  Is  correct  for  the  header 
of  a  MOSAIQUE  data  file.  This  header  Is  reproduced  as  the  header 
of  the  resulting  JOBNAME.ATI  file.  Thus,  It  will  be  the  header  for 
the  ATILA  Job. 


6.2.2.  THE  ENTRIES. 

The  ATILA  entries  ANALYSIS,  CLASS,  NLOAD,  PRECISION,  LCPDDC, 

MATERIALS,  GEOMETRY,  GEOMETRY  POLARIZATION,  REDUCTIONS, 
FREQUENCIES,  RADIATION  and  PRINTING,  If  provided  In  the  MOSAIQUE 
data  file,  are  simply  read  and  copied  In  the  resulting  JOBNAME.ATI 
file.  Obviously,  they  are  not  used  by  the  mesh  generation 

procedure.  The  END  entry  has  to  be  provided  to  MOSAIQUE.  As  well 
as  for  ATILA,  it  defines  the  end  of  the  entry  list.  It  will  also 
be  copied  in  the  resulting  JOBNAME.ATI  file. 

The  NODES  entry  Is  used  to  provide  the  coordinates  of  the  nodes 
which  define  the  super  elements.  Its  syntax  Is  exactly  the  same  as 
for  the  ATILA  data  file  and  the  node  coordinates  have  to  be  listed 
in  the  node  number  ascending  order  associated  to  the  gross  mesh. 
Afte”  the  MOSAIQUE  run,  this  NODES  entry  will  be  copied  In  the 

JOBNAME.ATI  file  with  the  coordinate  list  which  is  associated  with 
the  generated  nodes.  The  new  node  numbering  will  be  the  node 
numbering  provided  by  this  list.  It  has  to  be  noted  that  the 
numbering  of  the  nodes  belonging  to  the  Initial  gross  mesh  Is 
modified  during  the  generation.  Moreover,  as  well  as  for  ATILA, 
nodes  have  to  be  set  in  the  compulsory  order:  solid,  fluid, 

radiating.  Finally,  the  entry  NEWAXES  and  SCALE  can  be  used  to 
define  the  node  coordinates  for  the  gross  mesh  as  well  as  they  are 
used  for  ATILA. 

As  described  in  chapter  3,  the  ELEMENTS  entry  is  divided  In 
sets,  each  set  being  related  to  a  given  element  type  .  Successive 
sets  have  to  be  separated  by  a  dummy  line  and  the  whole  entry  must 
also  be  closed  by  a  dummy  line.  However,  with  respect  to  ATILA, 
the  content  of  a  set  l»  modified  in  MOSAIQUE  by  replacing  each 
topology  line  of  the  type: 

Nil  N12  N13  .  .NIP 

by  a  set  of  lines  of  the  type: 

Nil  N12  N13  .  .NIP 

-1  ITYP  NK  NL  IN  NN  I 

1X1  1X2  1X3  ...  .  IXNK 

1 Y 1  IY2  IY3  .  .  .  I YNL 
IZ1  IZ2  IZ3 . IZMi 


.‘-V-r. 


where : 


-  Nil,  N12,  N13,  ....  NIP  are  the  actual  numbers  of  the 

corresponding  super  element  nodes, 

-  the  second  line  (-1  ITYP  NK  NL  NM  NN  I)  provides  compulsory 
Information  about  the  splitting  of  this  super  element, 

-  ITYP  refers  to  a  given  frame  (see  section  6.3), 

-  NK  Is  the  number  of  generated  elements  along  the  OX  natural  axis 
of  the  super  element, 

-  NL  Is  the  number  of  generated  elements  along  the  OY  natural  axis 
of  the  super  element  (If  necessary), 

-  NM  Is  the  number  of  generated  elements  along  the  OZ  natural  axis 
of  the  super  element  (if  necessary), 

-  NN  Is  a  dummy  field,  kept  for  special  uses, 

-I  Is  a  control  Integer  which  must  be  equal  to  0  If  the 
splittings  are  regular  and  1  if  at  least  one  splitting  Is 
irregular , 

-  1X1,  1X2,  1X3,  ....  I XNK  define  the  splitting  along  the  OX 
natural  axis  of  the  super  element. 

-  IY1 ,  IY2,  IY3 .  I YNL  and  1 21 ,  IZ2,  IZ3 .  IZNM  provide 

the  same  Information,  respectively  for  the  OY  and  OZ  natural 
axes  of  the  super  element. 


Then,  along  the  OX  natural  axis  of  the  super  element,  the 
corresponding  side  Is  divided  Into  NK  different  parts,  the  lengths 
of  which  are  respectively  equal  to: 

IX1.L/IX,  IX2.L/IX,  1X3 .L/ IX . IXNK.L/IX 

where  L  Is  the  total  length,  IX  Is  given  by: 

IX  -  1X1  +  1X2  +  1X3  +  ...  +  I XNK 

and  1X1,  1X2,  1X3,  . ...IXNK  have  to  be  Integers.  The  same 
splitting  rule  Is  applied  for  the  OY  and  OZ  natural  directions.  If 
1  Is  equal  to  0,  the  three  lines  which  are  devoted  to  the 
description  of  the  splitting  are  not  needed.  If  I  Is  equal  to  1, 
these  three  lines  must  be  provided,  even  if  one  or  two  of  them 
correspond  to  a  regular  splitting. 

In  the  super  element  topology,  a  mid  side  node  can  always  be 
omitted  If  the  corresponding  side  Is  straight.  Nevertheless,  If  a 
topology  line  Is  Incomplete,  It  has  to  be  completed  with  zeros,  at 
least  until  the  last  defined  node.  As  an  example,  If  the  mid  side 
nodes  12,  14  and  21  can  be  omitted  In  the  topology  line  1  3  5  9  12 
14  16  21,  this  line  can  be  written: 

125900  16  0  or  125900  16 

In  opposition  to  the  ATILA  case,  when  solid  and  fluid  elements  are 
used  together,  the  elements  must  be  provided  In  the  order:  solid. 
Interface,  fluid,  damping.  This  order  can  be  modified  by  the  user 
In  the  J0BNAME .AT I  file  after  the  M0SAI0UE  run,  if  needed. 
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Finally,  two  new  entries  are  available  for  MOSAIQUE  with 
respect  to  ATILA.  The  first  Is  used  to  provide  Information  related 
to  the  generation  of  the  damping  elements  and  must  be  written  in 
the  following  form : 

SORTFLU ID 
NFLU  NAC  NSH 

where  NFLU  Is  the  actual  number  of  the  first  node  belonging  to  the 
boundary  surface  of  the  gross  mesh,  NAC  Is  the  actual  number  of 
the  node  which  Is  the  boundary  surface  centre  (acoustic  centre) 
and  NSH  Is  a  dummy  field  kept  for  special  uses.  If  the  boundary 
surface  centre  does  not  belong  to  the  mesh  and  cannot  be 
identified  with  a  node,  NCA  has  to  be  equal  to  zero  and  the  origin 
of  the  global  axes  must  be  this  centre.  This  entry  Is  compulsory 
if  damping  elements  are  used.  The  second  new  entry  Is  used  to 
modify  the  size  of  the  2D  or  3D  window  which  allows  the 
elimination  of  redundant  nodes.  It  has  to  be  used  If  redundant 
nodes  are  still  present  In  the  JOBNAME.ATI  file  or  If  actual  nodes 
have  been  deleted.  It  has  to  be  written  under  the  form.- 

MERGE 

XM 

where  XM  will  be  the  new  size  of  the  window  (XM  is  a  length  and 
must  be  provided  In  meters).  This  entry  has  only  to  be  used  If  the 
above  mentioned  problems  appear,  which  could  be  the  case  if  the 
mesh  Is  very  Irregular.  Then  the  best  way  of  choslng  the  XM  value 
Is  a  trial  and  error  method. 


In  fact,  the  element  type  Is  related  simultaneously  to  the 
super  element  and  the  generated  finite  elements.  This  situation  Is 
simple  when  the  super  element  and  the  generated  elements 
correspond  to  the  same  geometrical  shape  (splitting  of  a 
quadrilateral  super  element  Into  QUAD08E  finite  elements, 
splitting  of  a  triangular  super  element  Into  AXIS06P  finite 
elements  ...)  but  can  lead  to  severe  difficulties  If  geometrical 
shapes  are  different  (splitting  of  a  quadrilateral  super  element 
into  TRIA06E  finite  elements,  as  an  example).  In  this  lest  case, 
the  user  must  provide  an  element  type  which  corresponds  tc  the 
geometrical  shape  of  the  super  element.  Then,  MOSAIQUE  writes  down 
this  wrong  element  type  In  the  JOBNAME.ATI  file  and  the  user  has 
to  modify  this  type  before  using  the  data  file.  Work  Is  in 
progress  to  avoid  this  problem  In  the  next  ATILA  version. 


6.2.3.  THE  DATA. 

Data  can  be  simply  added  to  the  entries  In  the  JOBNAME.ATI 
file,  after  the  automatic  mesh  generation.  This  procedure  Is 
classical,  since  the  node  numbers  which  are  necessary  to  write 
down  the  data  are  only  available  after  the  generation.  However,  In 
special  cases  (multiple  runs  with  the  same  topology  and  different 
dimensions,  for  example),  the  user  may  wish  to  write  down  these 
data  in  the  JOBNAME.MOS  file.  Then,  the  rules  are  the  same  as  for 
an  ATILA  data  file  and  the  corresponding  lines  are  simply  read  by 
MOSAIQUE  and  copied  Into  the  JOBNAME.ATI  file. 
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6.3  THE  GENERATING  FRAMES 


The  splitting  of  a  super  element  In  simple  finite  elements  can 
be  realised  in  various  ways,  depending  upon  the  shape  of  this 
super  element  and  the  type  of  the  requested  finite  elements 
(quadrilateral  2D  elements  with  8  nodes,  20  node  hexaedras,  6  node 
ax Isymmetr ica I  Interfaces  ...).  The  available  types  of  splitting, 
named  frames  in  the  previous  sections,  are  marked  by  the  ITYP 
control  integer  and  described  In  the  following  lines. 


ITYP  -  1 

A  linear  super  element  Is  split  Into  2  node  linear  finite  elements 
(SPRI02E)  following  the  diagram: 


.0...0...0...2 

NK  elements  (NK  -  2) 


The  super  element  Is  defined  by  nodes  1  and  2  and  the  control 
integers  NL  and  NM  are  dummy. 


ITYP  -  2 

A  linear  super  element  Is  split  Into  3  node  damping  elements 
(AXIS03M  or  AXIS03D)  following  the  diagram: 


1...0...0...0...2 


NK  elements  (MC  -  2) 


The  super  element  Is  defined  at  least  by  the  nodes  1  and  2  but 
generally  also  by  its  mid  node,  due  to  the  curved  shape  of  the 
damping  boundary.  The  control  Integers  NL  and  NM  are  dummy. 


ITYP  -  3 


A  linear  super  element  Is  split  Into  6  node  Interface  elements 
(AXIS06I)  following  the  diagram: 


Fluid  domain 

0...0...0...0...0 

1...0...0...0...2 

Solid  doma I n 

NK  elements  (NK  -  2) 


The  super  element  Is  defined  by  the  nodes  1  and  2.  Its  mid  node  Is 
required  only  If  the  Interface  Is  curved.  The  control  Integers  NL 
and  NM  are  dummy.  The  orientation  of  the  normal  to  the  Interface 
from  the  fluid  to  the  solid  being  compulsory  for  the  generated 
Interface  finite  elements,  the  orientation  of  the  super  element  Is 
also  compulsory.  Nodes  1  and  2  are  defined  In  the  above  order  If 
the  Oz  natural  axis  Is  the  outward  normal  to  this  paper  sheet. 
Then,  the  Ox  natural  axis  is  oriented  toward  the  left  (attention 
must  be  devoted  to  this  point,  since  this,  numbering  rule  which 
applies  to  the  super  element  Is  the  opposite  of  the  rule  which 
applies  to  the  finite  elements). 


ITYP  -  4 

A  quadrilateral  super  element  Is  split  Into  8  node  quadrilateral 
elements  ( QUAD08E ,  PLAT08E ,  FACE08E,  AXIS08P,  QUAD08P,  AXIS08F  or 
QUAD08M)  following  the  diagram: 


NL 

elements 
(NL  -  2) 


NK  elements  (NK  -  3) 
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The  super  element  Is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  If  the  sides  are  curved.  The  control 
Integer  NM  is  dummy. 


ITYP  -  5 

A  quadrilateral  super  element  Is  split  Into  6  node  triangular 
elements  (TRIA06E.  PLAT06E,  AXIS06P,  TRIA06P,  AXIS06F  or  TRIA06M) 
following  the  diagram-. 


3...  0  .  «  .  0  .  *  •  0  »  •  .  4 


0  0 


0  0 


NL 

elements 
(NL  -  3) 


0..0..0..0..0..0 


0  0 


0  0 


.  0  ...0...  0  ...0 


0  0 


0  0 


0  0 


NK  elements  (NK  -  3) 


The  super  element  Is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  If  the  sides  are  curved.  The  control 
integer  NM  is  dummy.  This  case  Is  related  to  the  comments  at  the 
end  of  section  6.2.2  and  care  must  be  paid  to  the  element  type  (as 
an  example,  the  element  type  can  be  QUAD08E  In  the  MOSAIQUE  data 
file  and  TRIA06E  In  the  ATILA  data  file). 


ITYP  -  6 

A  triangular  super  element  Is  split  Into  6  node  triangular 
elements  (TRIA06E,  PLAT06E ,  AXIS06P,  TRIA06P,  AXIS06F  or  TRIA06M) 
following  the  diagram: 


v.v  v.v 
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The  super  element  Is  defined  by  the  nodes  1,  2  and  3.  M 
nodes  are  only  required  if  the  sides  are  curved.  The 
Integer  NM  is  dummy. 


ITYP  -  9 

A  quadrilateral  super  element  Is  split  Into  16  node  quadr I 
interface  elements  (QUAD16I)  following  the  diagram: 


NK  elements  (NK  -  2) 


id  side 
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The  super  element  Is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  If  the  sides  are  curved.  The  control 
Integer  NM  Is  dummy.  Attention  must  be  paid  to  the  orientation  of 
this  super  element,  the  normal  to  the  Interface  being  always 
oriented  from  the  fluid  toward  the  solid. 


ITYP  -  10 

An  hexaedric  super  element  Is  split  into  20  node  hexaedrlc  finite 
elements  (HEXA20E,  HEXA20P  or  HEXA20F)  following  the  diagram: 
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The  super  element  Is  defined  by  the  nodes  1,  2,  3,  4,  5,  6,  7  and 
8.  Mid  side  nodes  are  only  required  If  the  sides  are  curved. 
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ITYP  -  11 


A  prismatic  super  element  with  a  triangular  basis  Is  split  Into  15 
nodes  prismatic  finite  elements  (PRIS15E,  PRIS15P  or  PRIS15F) 
following  the  diagram: 
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The  super  element  Is  defined  by  nodes  1,  2,  3,  4,  5  and  6.  Mid 
side  nodes  are  only  required  If  the  sides  are  curved. 


ITYP  -  12 


An  hexaedrlc  super  element  Is  split  Into  15  node  prismatic  finite 
elements  (PRIS15E,  PRIS15P  or  PRIS15F)  following  the  diagram: 
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The  super  element  Is  defined  by  the  nodes  1,  2,  3,  4,  5,  6,  7  and 
8.  Mid  side  nodes  are  only  required  If  the  sides  are  curved.  This 
case  Is  related  to  the  comments  at  the  end  of  section  6.2.2  and 
care  must  be  paid  to  the  element  type  (as  an  example,  the  element 
type  can  be  HEXA20E  In  the  MOSAIQUE  data  file  and  PRIA15E  In  the 
ATILA  data  file). 


I TYP  -  14 

A  linear  super  element  Is  split  into  3  node  ax  I symmetr I ca I  shell 
finite  elements  (SHEL03E).  The  Information  related  to  ITYP  -  2  Is 
also  correct  for  this  case  but  the  numbering  of  the  generated 
nodes  and  the  topology  of  the  generated  elements  take  account  of 
the  specific  numbering  given  for  this  shell  In  section  4.3 
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